AD-A239  018 


July  15,  1991 


N00140- 88-0677 


Investigation  of  Hot  Streak  Migration  and 
Film  Cooling  Effects  on  the  Heat  Transfer 
in  Rotor/ Stator  Interacting  Flows 


Quarterly  Technical  Progress  Report  For 
April  1  through  June  30,  1991 

Prepared  for  NAVAIR 

by 

Daniel  J.  Dorney  and  Roger  L.  Davis 
Theoretical  &  Computational  Fluid  Dynamics 
United  Technologies  Research  Center 

Introduction 

The  specific  objectives  of  this  effort  are  to  extend  the  NASA-Ames  Rai  [1,  2]  three- 
dimensional  rotor /stator  interaction  code,  R0T0R3,  to  investigate  hot  streak  migration 
and  film  cooling  effects  on  the  passage  flow  and  blade  surface  heat  transfer  for  an  axial  flow 
turbine  stage.  These  objectives  are  part  of  an  overall  plan  to  extend  the  capabilities  of  this 
numerical  procedure  and  to  determine  the  potential  of  this  technique  to  impact  the  design 
of  future  rotating  turbomachinery  components. 


The  principal  benefits  that  will  result  from  this  effort  are: 

•  A  diagnostic  analysis  and  code  which  can  be  used  to  help  design  turbine  blades  with 
increased  efficiency  and  reduced  cooling  requirements. 

•  An  open  literature  demonstration  on  the  use  of  Computational  Simulation  and  Scien¬ 
tific  Visualization  for  gaining  insight  into  complex  turbomachinery  flows. 

•  Acceleration  of  the  transition  of  large-scale  computational  analyses  to  the  turbomachin¬ 
ery  design  process  through  joint  involvement  between  UTRC  and  Pratt  and  Whitney 
turbine  engineers. 

•  A  useful  and  proven  analytical  tool  that  can  be  used  with  existing  and  future  engine 
design  procedures. 


The  specific  tasks  associated  with  this  effort  are: 
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•  The  R0T0R3  analysis  will  be  extended  by  incorporating  a  three-dimensional  hot 
streak  at  the  stator  inlet  to  numerically  simulate  the  flow  which  exits  a  combustor. 
The  effects  of  the  combustor  hot  streak  on  turbine  stator  and  rotor  passage  flows  and 
blade  surface  heat  transfer  characteristics  will  be  investigated.  In  addition,  the  RO¬ 
TORS  code  will  be  extended  to  include  a  film  cooling  model  in  an  effort  to  understand 
the  effects  of  film  cooling  on  stator  and  rotor  blade  surface  heat  transfer  and  on  the 
flow  through  the  downstream  rotor.  Both  the  hot  streak  and  heat  transfer  investiga¬ 
tions  will  be  validated  against  data  obtained  from  the  UTRC  Large  Scale  Rotating  Rig 
(LSRR)  experiments. 

•  The  UTRC  Scientific  Visualization  Laboratory  will  provide  a  means  of  exhibiting  the 
important  flow  features.  Advanced  graphical  representations  are  essential  for  reduc¬ 
ing  the  massive  amounts  of  numerical  data  that  result  from  a  calculation  of  three- 
dimensional  viscous  flows  through  multi-blade  row  turbomachinery  passages.  Scientific 
Visualization  will  be  used  to  display  details  of  the  hot  streak  and  film  cooling  flow  paths 
as  well  as  to  provide  information  on  blade  surface  heat  transfer  and  mainstream  flow 
temperature.  The  UTRC  Scientific  Visualization  Lab  will  provide  animated  recordings 
for  use  by  the  Navy  and  others. 


Quarterly  Technical  Progress  Summary 

Task  I:  Modify  Code  Boundary  Conditions 

No  work  was  planned  for  Task  I  during  this  quarter 


Task  II:  Verify  and  Improve  Code 

Two  technical  papers  were  presented  at  meetings  during  the  past  quarter  which  encom¬ 
pass  work  completed  under  Tasks  II  and  III.  The  first  technical  paper,  entitled  “Navier- 
Stokes  Analysis  of  Turbine  Blade  Heat  Transfer  and  Performance,”  was  presented  at  the 
AGARD  Symposium  on  CFD  Techniques  for  Propulsion  Applications  in  San  Antonio,  Texas 
during  the  week  of  May  27.  The  objective  of  this  paper  was  to  establish  the  grid  density 
requirements  for  accurately  predicting  blade  surface  and  endwall  heat  transfer  in  a  turbine 
passage.  Simulations  of  transitional  flow  through  the  Langston  cascade  [3,  4]  with  a  thin  inlet 
boundary  layer  and  heat  transfer  were  performed.  It  was  determined  that  for  this  particular 
geometry  approximately  900,000  computational  grid  points  are  required  to  produce  nearly 
grid-independent  flow  solutions.  In  addition,  the  study  showed  that  overall  aerodynamic  tur- 
bomachinery  blade  row  performance  can  be  predicted  within  2%  using  fine  computational 
grids  and  the  Baldwin-Lomax  algebraic  turbulence  model.  The  second  technical  paper,  en¬ 
titled  “Two-Dimensional  Inlet  Temperature  Profile  Attenuation  in  a  Turbine  Stage,”  was 
presented  at  the  36th  ASME  International  Gas  Turbine  and  Aeroengine  Congress  and  Ex - 
position  in  Orlando,  Florida,  during  the  week  of  June  3.  The  purpose  of  this  paper,  which 
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was  prepared  under  a  parallel  Corporate  Sponsored  Research  Program,  was  to  determine  if 
a  two-dimensional  Navier-Stokes  procedure  could  be  used  to  predict  unsteady  viscous  flow 
through  a  1-stator/ 1-rotor  configuration  with  a  planar  hot  streak  introduced  at  the  stator 
inlet.  Using  a  new  method  for  correlating  the  predicted  results  with  the  experimental  data, 
the  computational  procedure  accurately  predicted  the  rotor  airfoil  pressure  surface  heating 
due  to  the  hot  streak.  This  investigation  also  demonstrated  the  importance  of  hot  streak 
shape  on  the  predicted  results.  In  this  investigation,  the  predicted  results  obtained  with  a 
two-dimensional  Navier-Stokes  procedure  were  compared  with  experimental  data  for  a  planar 
(i.e.  two-dimensional)  hot  streak,  whereas  in  previous  two-dimensional  numerical  hot  streak 
simulations  the  computational  results  were  compared  with  experimental  data  obtained  for  a 
circular  (i.e.  fully  three-dimensional)  hot  streak.  Thus,  for  certain  hot  streak  shapes,  two- 
dimensional  numerical  simulations  can  be  used  in  place  of  more  expensive  three-dimensional 
numerical  simulations,  giving  designers  an  accurate  and  efficient  tool  for  optimizing  cool¬ 
ing  systems.  The  papers  are  attached  to  this  progress  report  as  they  are  related  to  work 
completed  to  date  under  Tasks  II  and  III. 

Task  III:  Perform  Numerical  Experiments 

The  following  numerical  experiment  is  being  performed  using  the  ROTOR3  code: 

•  A  1-stator/l-rotor  hot  streak  simulation  with  film  cooling  and  heat  transfer  has  been 
partially  completed.  In  this  investigation,  the  number  and  location  of  the  film  cooling 
slots  has  been  varied  to  obtain  both  the  qualitative  and  quantitative  characteristics  of 
the  effectiveness  of  different  film  cooling  schemes. 

The  computational  grid  topology  used  in  the  three-dimensional  numerical  simulation  is 
similar  to  those  used  in  a  previous  simulation  [5].  The  stator  grid  system  is  constructed  with 
101  x  21  grid  points  in  each  spanwise  O-grid  and  58  x  31  grid  points  in  each  spanwise  H-grid. 
The  rotor  grid  system  is  constructed  with  101  x  21  grid  points  in  each  spanwise  O-grid  and 
60  x  31  grid  points  in  each  spanwise  H-grid.  A  total  of  51  O-H  grid  planes  are  distributed 
in  the  spanwise  direction.  The  rotor  region  has  a  tip  clearance  grid  system  that  contains 
101  x  11  grid  points  in  each  of  7  spanwise  locations.  A  total  of  410,677  grid  points  are  used 
in  the  three-dimensional  simulation.  A  wall  spacing  of  1.27  x  10~5  meters  (y+  w  4)  is  used 
in  the  blade-to-blade  direction,  while  a  wall  spacing  of  7.62  x  10-4  meters  ( y+  ss  18)  is  used 
in  the  spanwise  direction. 

In  the  first  portion  of  the  simulation,  a  heat  flux  was  applied  to  the  rotor  surface  with¬ 
out  film  cooling.  The  heat  flux  was  specified  to  yield  a  surface  temperature  approximately 
3-4  degrees  less  than  the  local  adiabatic  surface  temperature.  Figure  1  illustrates  the  ro¬ 
tor  surface  mid-span  time-averaged  temperature  coefficient  distributions  for  the  adiabatic 
simulation  [5]  and  the  simulation  with  rotor  surface  heat  transfer,  where  the  time-averaged 


3 


Figure  1:  Rotor  surface  mid-span  time-averaged  temperature  distributions  with  and  without 
heat  transfer 


temperature  coefficient  is  defined  as: 

T  —  T 

Ct  =  t - ~ZPt  (!) 

1  av9rle  1 1 

and  T\  is  the  free  stream  temperature  outside  the  influence  of  the  hot  streak,  T  is  the  local 
time-averaged  surface  temperature,  and  Tavgrlt  is  the  time-averaged  surface  temperature  at 
the  rotor  mid-span  leading  edge.  The  addition  of  a  heat  flux  is  seen  to  have  little  effect  on 
the  time-averaged  rotor  surface  temperature  distribution.  This  is  expected  considering  the 
relatively  small  specified  heat  flux  and  the  high  local  time-averaged  surface  temperatures 
due  to  the  hot  streak.  A  more  negative  heat  flux  value  may  reduce  the  influence  of  the  hot 
streak  on  the  surface  temperature,  but  may  also  have  adverse  effects  on  the  stability  of  the 
computational  procedure.  Figure  2  shows  time-averaged  temperature  contours  for  the  rotor 
pressure  surface  with  heat  transfer.  The  time-averaged  temperature  contours  for  this  case 
are  nearly  identical  to  those  obtained  from  the  adiabatic  hot  streak  simulation. 

In  the  second  portion  of  the  simulation,  fluid  was  injected  through  holes  in  the  rotor 
pressure  surface  at  a  velocity  equal  to  0.2  times  that  of  the  free  stream  velocity  and  at  an 
angle  of  25  degrees  with  respect  to  the  local  axial  direction  surface  tangent.  The  temperature 
of  the  injected  fluid  was  specified  to  be  520°  Rankine  (2%  lower  than  the  free  stream  and  22% 
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TIME-AVERAGED  TEMPERATURE  CXMOURS 
PRESSURE  SURFACE 


LEADING  EDGE 


Figure  2:  Time-averaged  temperature  contours  for  rotor  pressure  surface  with  heat  transfer 
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Figure  3:  Rotor  surface  mid-span  time-averaged  temperature  distributions  for  the  adiabatic, 
heat  transfer  only,  and  two  row  film  cooling/heat  transfer  simulations 


lower  than  the  hot  streak)  and  the  density  ratio  was  chosen  such  that  the  static  pressures 
at  the  injection  locations  were  equal  to  the  time-averaged  static  pressures  in  the  absence  of 
film  cooling.  The  injection  holes  (ports)  were  located  at  20%  and  40%  of  the  axial  chord  and 
extended  from  20%  of  the  span  to  80%  of  the  span  at  approximately  5%  spanwise  intervals. 
Similar  to  previous  two-dimensional  simulations  [6],  two-adjacent  grid  points  were  used  to 
model  an  injection  hole.  A  total  of  22  injections  holes  were  used  in  the  first  portion  of 
the  simulation.  The  surface  heat  flux  was  again  specified  to  yield  a  surface  temperature 
approximately  3-4  degrees  less  than  the  local  adiabatic  surface  temperature. 

Figure  3  illustrates  the  rotor  surface  mid- span  time- averaged  temperature  coefficient 
distributions  for  the  adiabatic  simulation,  the  simulation  with  rotor  surface  heat  transfer 
only,  and  the  simulation  with  two  film  cooling  rows  and  heat  transfer.  A  substantial 
reduction  of  the  time-averaged  temperature  on  the  rotor  pressure  surface  is  observed  in 
the  vicinity  of  the  film  cooling  injection  holes.  As  the  cooling  fluid  convects  downstream 
and  away  from  the  airfoil  surface,  the  time-averaged  surface  temperature  increases  due  to 
the  influence  of  the  hot  streak.  In  this  three-dimensional  simulation,  the  time-averaged 
surface  temperature  in  the  region  between  film  cooling  injection  holes  and  downstream  of 
the  last  film  cooling  hole  is  seen  to  increase  much  more  rapidly  than  that  observed  in  similar 
two-dimensional  simulations  [6].  This  may  be  partially  due  to  the  interaction  between  the 
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secondary  and  wall  layer  flows  in  the  rotor  passage,  which  cause  the  cooler  injected  fluid 
to  mix  very  rapidly  with  the  hot  streak  fluid.  Another  cause  of  the  temperature  increase 
may  be  that  the  current  injection  angle  (25  degrees)  is  greater  than  that  used  in  the  two- 
dimensional  simulations  (10  degrees).  While  it  appears  to  be  beneficial  to  reduce  the  injection 
angle,  realistically  structural  and  manufacturing  considerations  would  limit  the  actual  range 
of  possible  injection  angles.  The  current  injection  angle  was  chosen  to  be  commensurate 
with  those  found  in  actual  gas  turbines.  In  spite  of  placing  film  cooling  injection  holes  at 
20%  and  40%  axial  chord  locations,  the  rapid  mixing  between  the  cooling  fluid  and  the  hot 
streak  causes  an  increase  in  the  temperature  of  the  pressure  surface  near  the  trailing  edge 
region.  The  time-averaged  temperature  on  the  rotor  suction  surface  is  unaffected  since  the 
film  cooling  holes  were  placed  only  on  the  pressure  surface. 

Figure  4  shows  time-averaged  temperature  contours  for  the  rotor  pressure  surface  with 
2  rows  of  film  cooling  holes  and  heat  transfer.  The  film  cooling  locations  are  characterized 
by  the  dark  blue  (low  temperature)  contours.  While  discrete  film  cooling  jets  can  be  seen 
inboard  of  approximately  35%  span  and  outboard  of  about  65%  span,  the  cooling  fluid  near 
mid-span  appears  to  be  contained  in  one  continuous  jet.  As  noted  above,  the  film  cooling 
injection  holes  were  spaced  approximately  every  5%  of  the  span.  In  the  mid-span  region, 
the  computational  grid  lines  were  also  spaced  approximately  5%  of  the  span  apart,  resulting 
in  6  film  cooling  holes  being  placed  at  six  consecutive  grid  points.  Figure  4  shows  that  the 
combination  of  cooling  fluid  injection  angle  and  secondary /boundary  layer  flow  causes  rapid 
mixing  of  the  cooling  fluid  and  allows  the  hot  streak  to  increase  the  rotor  pressure  surface 
temperature  aft  of  mid-chord  region. 

In  an  effort  to  eliminate  the  high  time-averaged  temperatures  aft  of  mid-chord,  a  third 
row  of  film  cooling  holes  was  added  at  approximately  60%  axial  chord.  Thus,  in  this  portion 
of  the  simulation  a  total  of  33  film  cooling  holes  (each  consisting  of  two  adjacent  grid  points) 
were  used.  In  addition,  the  injection  velocity  was  increased  to  0.4  times  that  of  the  free 
stream  velocity.  Figure  5  illustrates  the  rotor  surface  mid-span  time-averaged  temperature 
coefficient  distributions  for  the  adiabatic  simulation,  the  simulation  with  rotor  surface  heat 
transfer  only,  and  the  simulation  with  three  film  cooling  rows  and  heat  transfer.  In  this 
simulation,  the  combination  of  a  greater  injection  velocity  and  a  third  row  of  film  cooling 
injection  holes  has  effectively  eliminated  most  of  the  high  temperature  regions  on  the  rotor 
pressure  surface.  The  high  temperature  region  aft  of  mid-chord  has  been  reduced,  as  has 
the  high  temperature  regions  between  the  film  cooling  injection  holes.  A  small  spike  in 
the  time-averaged  temperature  distribution  is  noticeable  just  upstream  of  the  first  row  of 
film  cooling  holes  and  just  downstream  of  the  first  and  second  rows  of  film  cooling  holes. 
These  spikes  may  be  due  to  the  sudden  change  from  a  no-slip  to  a  transpiration  boundary 
condition,  resulting  in  sudden  changes  in  the  density  and  velocity.  The  magnitude  of  these 
spikes  increase  as  the  film  cooling  velocity  is  increased,  thus  creating  larger  gradients  in 
the  flow  variables.  One  possible  technique  for  avoiding  spikes  is  to  gradually  increase  and 
decrease  the  film  cooling  parameters  of  adjacent  computational  nodes. 
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TIME-AVERAGED  TEMPERATURE  CONTOURS 
PRESSURE  SURFACE 


Figure  4:  Time-averaged  temperature  contours  for  rotor  pressure  surface  with  two  rows  of 
film  cooling  holes  and  heat  transfer 
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Figure  5:  Rotor  surface  mid-span  time-averaged  temperature  distributions  for  the  adiabatic, 
heat  transfer  only,  and  3  row  film  cooling/heat  transfer  simulations 
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Figure  6  shows  time-averaged  temperature  contours  for  the  rotor  pressure  surface  with  3 
rows  of  film  cooling  holes  and  heat  transfer.  The  film  cooling  locations  are  again  characterized 
by  the  dark  blue  (low  temperature)  contours.  Comparing  Figs.  4  and  6  shows  that  increasing 
the  injection  velocity  extends  the  range  of  effectiveness  of  the  cooling  fluid.  Figure  6  also 
illustrates  that  the  addition  of  a  third  row  of  film  cooling  holes  has  significantly  decreased 
the  time-averaged  rotor  pressure  surface  temperature  aft  of  the  mid-chord  point.  A  fourth 
row  of  film  cooling  holes  placed  at  approximately  80%  chord  would  probably  eliminate  the 
lingering  remnants  of  the  hot  streak  in  the  trailing  edge  region. 

Task  IV:  Utilize  Scientific  Visualization 

A  video  animation  of  the  predicted  results  from  a  two-dimensional  3-stator/4-rotor/l- 
hot  streak  simulation  with  film  cooling  and  heat  transfer  was  created.  The  video,  which 
illustrated  static  temperature  contours,  was  constructed  using  80  time  segments  obtained 
during  one  passing  cycle  of  the  rotor.  The  video  was  used  to  study  the  migration  of  the 
cooling  fluid,  which  was  injected  from  the  trailing  edge  of  the  stators  and  from  the  pressure 
surface  of  the  rotors. 


Work  Planned 

Task  I:  Modify  Code  Boundary  Conditions 
No  work  is  planned  for  Task  I. 

Task  II:  Verify  and  Improve  Code 
No  work  is  planned  for  Task  II. 

Task  III:  Perform  Numerical  Experiments 

The  three-dimensional  1-stator/Trotor/l-hot  streak  simulations  with  film  cooling  and 
heat  transfer  will  be  completed  and  the  predicted  unsteady  and  time-averaged  flow  fields 
will  be  analyzed  in  greater  detail. 

Task  TV:  Utilize  Scientific  Visualization 

Work  will  be  initiated  on  visualizing  and  animating  the  predicted  flow  field  for  the  three- 
dimensional  hot  streak  simulation  including  the  effects  of  film  cooling  and  heat  transfer. 
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TIME-AVERAGED  TEMPERATURE  CONTOURS 
PRESSURE  SURFACE 


Figure  6:  Time-averaged  temperature  contours  for  rotor  pressure  surface  with  three  rows  of 
film  cooling  holes  and  heat  transfer 
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Task  IV:  Reporting 


Work  on  the  final  technical  report  will  be  initiated. 
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ABSTRACT 

Experimental  evidence  has  shown  that  hot  streaks  which  are 
emitted  from  the  combustors  of  gas  turbines  are  often  largely 
responsible  for  the  burning  of  first  stage  turbine  blades.  Design¬ 
ers  have  attempted  to  counteract  the  effects  of  these  hot  streaks 
through  the  use  of  complex  internal  and  film  cooling  schemes. 
Unfortunately,  due  to  the  lack  of  accurate  predictive  tools  which 
account  for  temperature  non- uniformities  in  the  gas  path,  as  well 
as  a  lack  of  detailed  understanding  of  the  physical  mechanisms 
which  control  the  migration  and  accumulation  of  the  hot  streak 
gases,  turbine  blade  “hot  spots”  still  occasionally  occur.  In  an 
effort  to  increase  understanding  of  the  interaction  mechanisms 
between  combustor  hot  streaks  and  turbine  blade  heat  transfer, 
a  numerical  investigation  has  been  conducted  to  determine  if  a 
two-dimensional  solution  procedure  can  accurately  predict  rotor 
airfoil  surface  heating  for  flows  which  include  planar  hot  streaks. 
A  two-dimensional  Navier-Stokes  analysis  is  used  to  predict  un¬ 
steady  viscous  flow  through  a  1-stator/l-rotor  configuration  with 
a  planar  hot  streak  introduced  at  the  stator  inlet.  Comparison 
of  the  predicted  results  with  a  new  experimental  data  set  demon¬ 
strates  that  the  two-dimensional  numerical  procedure  can  be  used 
to  accurately  predict  time-averaged  rotor  pressure  surface  tem¬ 
peratures  for  flows  which  include  planar  hot  streaks. 

Nomenclature 

a  -  Speed  of  sound 

e  -  Specific  energy 

et  -  Total  energy 

J  -  Jacobian  of  transformation 
Af  -  Mach  number 

P  -  Static  pressure 

Pt  -  Prandtl  number 

R  -  Universal  gas  constant 


Re  - 

Free  stream  inlet  reference  Reynolds  number 

T 

Static  temperature 

u,  v  - 

x,  y  components  of  velocity 

U 

Rotor  velocity 

0  - 

Rotor  inlet  relative  flow  angle 

K 

Thermal  conductivity 

n,\  - 

First  and  second  coefficients  of  viscosity 

P 

Density 

r 

Shear  stress 

Subscript* 

c 

•  Cold  flow 

HS 

-  Hot  streak 

i 

-  Inviscid 

L 

-  Laminar  quantity 

t 

•  Stagnation  quantity 

T 

-  Turbulent  quantity 

V 

-  Viscous 

w 

-  Wall  value 

*,V 

-  First  derivative  with  respect  to  x  or  y 

xx,  VV 

-  Second  derivative  with  respect  to  x.  y 

1 

-  Inlet  quantity 

2 

-  Exit  quantity 

INTRODUCTION 

In  recent  years,  design  engineers  have  achieved  substantial 
increases  in  turbine  stage  efficiencies  by  increasing  combu*tur 
operating  temperatures.  The  increased  temperature  of  the  cun' 
bustor  exit  flow  has  necessitated  the  implementation  of  complex 
cooling  schemes  to  avert  burning  of  first  stage  turbine  airfoil*  In 
an  effort  to  optimize  the  various  cooling  schemes,  engineer*  ha'* 
simulated  (both  experimentally  (Butler  ef  of.  1989)  and  ■  < -n . p o 


tationally  (Krouthen  and  Giles,  1988,  Dorney  et  al ,  1990,  Taka¬ 
hashi  and  Ni,  1990,  Sharma  et  al  1990,  Rai  and  Dring,  1990))  the 
“hot  streaks”  which  exit  the  combustor  and  migrate  through  the 
turbine  stage.  The  goal  of  these  simulations  has  been  to  isolate 
the  various  factors  associated  with  combustor  hot  streaks  which 
lead  to  excessive  temperatures  and  burning  of  first  stage  turbine 
blade  pressure  surfaces. 

The  work  of  Butler  et  al  (1989)  represents  a  detailed  experi¬ 
mental  investigation  of  hot  streak  migration.  In  this  experiment, 
a  three-dimensional  hot  streak  was  introduced  into  the  stator  in¬ 
let  of  a  turbine  stage  through  a  circular  pipe.  The  hot  streak  was 
seeded  with  CO2,  and  the  path  of  the  hot  streak  was  tracked  by 
sampling  the  CO2  concentrations  at  various  locations  in  the  tur¬ 
bine  stage.  According  to  Butler  et  al  (1989),  the  experimental 
simulations  indicate  that  hot  streak  migration  paths  and  rotor 
surface  heating  are  controlled  by  a)  the  secondary  flows  in  the 
rotor  passage,  and  b)  the  segregation  of  hot  and  cold  gases  caused 
by  the  difference  in  the  rotor  inlet  flow  angle  between  the  hot  and 
cold  gases. 

Attempts  to  reproduce  the  circular  hot  streak  experimental 
results  of  Butler  et  al  (1989)  using  two-dimensional  computa¬ 
tional  techniques  have  only  been  partially  successful  (Krouthen 
and  Giles,  1988,  Dorney  et  al  1990,  Rai  and  Dring,  1990).  The 
discrepancies  between  the  predicted  two-dimensional  results  and 
experimental  data  are  probably  due  to  a  combination  of  fac¬ 
tors,  including  a)  the  effects  of  secondary  and  endwall  flows,  b) 
the  three-dimensional  shape  of  this  particular  experimental  hot 
streak,  c)  the  assumed  correlation  between  COa  (used  in  the 
experiments  to  track  the  movement  of  the  hot  streak)  and  tem¬ 
perature,  d)  the  hot  streak  to  rotor  passage  count,  and  e)  the 
ratio  of  the  hot  streak  temperature  to  that  of  the  surrounding 
flow. 

Some  progress  has  been  made  in  isolating  the  role  of  three- 
dimensional  and  viscous  effects  on  hot  streak  migration  (Dorney 
et  al  1990,  Takahashi  and  Ni,  1990,  Sharma  et  al  1990)  and 
turbine  rotor  heat  transfer.  Both  a  three-dimensional  Navier- 
Stokes  (Dorney  et  al  1990)  and  an  Euler  simulation  (Takahashi 
and  Ni,  1990,  Sharma  et  al  1990)  in  which  viscous  effects  were 
modelled  have  demonstrated  better  agreement  with  the  circular 
hot  streak  experimental  data  (Butler  et  al  1989)  than  previous 
two-dimensional  simulations  for  the  time-averaged  temperature 
on  the  rotor  pressure  surface. 

The  goal  of  the  current  research  effort  is  to  determine  if  two- 
dimensional  unsteady  rotor/stator/hot  streak  simulations  can 
provide  reliable  warnings  as  to  when  rotor  surface  temperatures 
become  too  high  due  to  combustor  hot  streaks.  If  two-dimensional 
simulations  can  be  shown  to  be  useful  in  providing  this  guid¬ 
ance,  then  the  need  to  execute  costly  and  time-consuming  three- 
dimensional  simulations  can  be  reduced  or  eliminated.  Compu¬ 
tational  results  for  a  two-dimensional  hot  streak  simulation  are 
compared  with  new  experimental  data  for  a  two-dimensional  pla¬ 
nar  hot  streak  which  was  introduced  across  the  span  of  a  turbine 
stage  stator  inlet.  An  important  mechanism  which  has  been  iden¬ 
tified  as  a  controlling  factor  on  the  surface  temperature  of  turbine 


rotors  is  hot  streak  shape.  By  understanding  the  strengths  and 
limitations  of  two-dimensional  hot  streak  simulations,  useful  in-  * 
sights  can  be  obtained  which  will  improve  turbine  blade  designs.  f 

Results  of  a  1-stator/l-rotor  hot  streak  migration  simulation 
of  the  Large  Scale  Rotating  Rig  (LSRR)  (Butler  et  al,  1989)  are 
compared  with  experimental  data  recently  obtained  for  a  planar 
hot  streak.  The  close  agreement  between  the  numerical  results 
and  the  experimental  data  suggests  that  two-dimensional  numer¬ 
ical  analyses  can  be  used  to  accurately  predict  rotor  pressure 
surface  temperature  distributions. 


NUMERICAL  INTEGRATION  PROCEDURE 


The  governing  equations  of  motion  considered  in  this  study 
are  the  time  dependent  two-dimensional  Navier-Stokes  equations, 
which  can  be  written  in  Cartesian  coordinates  as: 
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For  the  present  application,  the  second  coefficient  of  viscosity  is 
calculated  using  Stokes'  hypothesis,  X  =  -2/3p.  The  equations 
of  motion  are  completed  by  the  perfect  gas  law. 

The  viscous  fluxes  can  be  simplified  by  incorporating  the  thin 
layer  assumption  (Baldwin  and  Lomax,  1978).  The  thin  layer 
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assumption  states  that  for  high  Reynolds  number  flows,  the  dif¬ 
fusion  terms  normal  to  a  solid  surface  will  be  much  greater  than 
*  those  parallel  to  the  surface.  In  the  current  study,  viscous  terms 
are  retained  in  the  direction  normal  to  the  blade  surface. 

To  extend  the  equations  of  motion  to  turbulent  flows,  an  eddy 
viscosity  formulation  is  used.  Thus,  the  effective  viscosity  and 
effective  thermal  conductivity  can  be  defined  as: 

H  =  HL+HT 


The  turbulent  viscosity,  /it,  is  calculated  using  the  two-layer 
Baldwin  and  Lomax  (1978)  algebraic  turbulence  model. 

The  numerical  procedure  for  the  two-dimensional  analysis 
consists  of  a  time  marching,  implicit,  third-order  spatially  ac¬ 
curate,  upwind,  finite  difference  scheme.  The  inviscid  fluxes 
are  discretized  according  to  the  scheme  developed  by  Osher 
(Chakravarthy  and  Osher,  1982).  The  viscous  fluxes  are  calcu¬ 
lated  using  standard  central  differences.  An  alternate  direction, 
approximate-factorization  technique  is  used  to  compute  the  time 
rate  changes  in  the  primary  variables.  In  addition,  an  inner  New¬ 
ton  iteration  is  used  to  increase  stability  and  reduce  linearization 
errors. 

GRID  GENERATION 

The  current  two-dimensional  Navier-Stokes  analysis  uses  mul¬ 
tiple  zonal  grids  to  discretize  the  rotor/stator  flow  field  and  facil¬ 
itate  relative  motion  of  the  rotor  (Rai  et  al  1988,1989).  A  com¬ 
bination  of  O-  and  H-grid  sections  are  generated  in  the  blade- to- 
blade  direction  extending  upstream  of  the  stator  leading  edge  to 
downstream  of  the  rotor  trailing  edge.  Algebraically  generated 
H-grids  are  used  in  the  regions  upstream  of  the  leading  edge, 
downstream  of  the  trailing  edge  and  in  the  inter-blade  region. 
The  O-grids,  which  are  body-fitted  to  the  surfaces  of  the  airfoils 
and  generated  using  an  elliptic  equation  solution  procedure,  are 
used  to  properly  resolve  the  viscous  flow  in  the  blade  passages  and 
to  easily  apply  the  algebraic  turbulence  model.  Computational 


H-GRID 


Figure  1:  2-D  computational  grid 


grid  lines  within  the  O-grids  are  stretched  in  the  blade-normal 
direction  with  a  fine  grid  spacing  at  the  wall.  Figure  1  illustrates 
the  grid  topology  used  in  the  current  simulations. 

BOUNDARY  CONDITIONS 

The  theory  of  characteristics  is  used  to  determine  the  bound¬ 
ary  conditions  at  the  stator  inlet  and  rotor  exit.  For  subsonic 
inlet  flow,  the  total  pressure,  tangential  velocity  component  and 
the  downstream  running  Riemann  invariant,  Rt  —  u  +  are 
specified  while  the  upstream  running  Riemann  invariant,  R2  = 
ti  -  is  extrapolated  from  the  interior  of  the  computational 
domain.  Inlet  flow  boundary  conditions  within  the  hot  streak 
region  are  updated  using  Riemann  invariants  where  the  velocity 
and  speed  of  sound  reflect  the  static  and  stagnation  temperature 
increase  in  the  hot  streak.  The  static  and  total  pressure  in  the  hot 
streak  are  assumed  to  be  equal  to  that  of  the  undisturbed  inlet 
flow,  corresponding  to  the  experimental  conditions.  For  subsonic 
outflow  the  pressure  is  specified  while  the  tangential  component 
of  velocity,  entropy,  and  the  downstream  running  Riemann  in¬ 
variant  are  extrapolated  from  the  interior  of  the  computational 
domain.  Absolute  no-slip  boundary  conditions  are  enforced  on 
stator  airfoil  surfaces,  while  relative  no-slip  boundary  conditions 
are  imposed  along  the  rotor  airfoil  surfaces.  Periodicity  is  en¬ 
forced  along  the  outer  boundaries  of  the  H-grids  in  the  circum¬ 
ferential  (y)  direction.  In  the  present,  study  the  flow  is  assumed 
to  be  adiabatic. 

Dirichlet  conditions,  in  which  the  time  rate  change  in  the 
vector  U  of  Eq.  (2)  is  set  to  zero,  are  imposed  at  the  patched 
boundaries  of  the  0-  and  H-grids  in  the  stator  and  rotor  regions. 
The  flow  variables  of  U  at  the  zonal  boundaries  are  explicitv 
updated  after  each  time  step  by  interpolating  values  from  the 
adjacent  grid.  Because  of  the  explicit  application  of  the  zonal 
boundary  conditions,  large  time  steps  necessitate  the  use  of  more 
than  one  Newton  iteration.  The  zonal  boundary  conditions  are 
non-conservative,  but  for  subsonic  flow  this  should  not  affect  the 
accuracy  of  the  final  flow  solution.  Further  information  describ¬ 
ing  the  implementation  of  the  boundary  conditions  can  be  found 
in  Rai  e<  of  (1985.1988,1989,1990). 

RESULTS 

A  numerical  investigation  of  hot  streak  migration  has  been 
conducted  in  which  the  predicted  results  from  a  two-dimensional 
Navier-Stokes  procedure  have  been  compared  with  experimental 
data.  The  turbine  airfoil  geometry  used  in  this  study  consisted 
of  the  first  stage  of  the  UTRC  Large  Scale  Rotating  Rig  ( LSRR  1 
(Dring  et  al ,  1982,  Butler  et  al  1989)  which  includes  22  stator 
airfoils  and  28  rotor  airfoils.  An  accurate  simulation  of  this  con 
figuration  would  require  at  least  11  stator  airfoils  and  14  rotor 
airfoils.  For  the  current  1-stator/l-rotor  simulation,  a  rescaling 
strategy  was  used  to  reduce  the  number  of  airfoils  to  1  stator 
and  1  rotor.  It  was  assumed  that  there  were  28  stator  airfoils 
and  28  rotor  airfoils  and  the  stator  was  scaled  down  by  the  lac 
tor  (22/28).  To  replicate  the  test  conditions  of  the  experimental 
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study,  a  65%  axial  gap  was  used  between  the  stator  and  rotor 
airfoils.  The  numerical  simulation  was  performed  at  the  exper¬ 
imental  stator  inlet  Mach  number  of  M  =  .051.  In  both  the 
experiment  and  the  numerical  simulation,  the  stator  vanes  were 
rotated  down  7  degrees  with  respect  to  the  tangential  direction 
relative  to  the  design  point  operating  conditions.  The  experimen¬ 
tal  rotor  rotation  speed  was  710  rpm.  The  Reynolds  number  in 
both  the  experiment  and  numerical  simulation  was  100,000  per 
inch  and  a  pressure  ratio  of  Pi/Pn  —  .9640  was  determined  from 
the  inlet  total  pressure  and  the  static  pressure  measured  in  the 
trailing-edge  plane. 

In  the  experimental  study,  one  hot  streak  was  introduced  in 
the  form  of  a  two-dimensional  jet  from  the  hub  to  the  tip  be¬ 
tween  two  stator  airfoils  of  the  LSRR  (i.e.  one  hot  streak  for  the 
entire  wheel).  The  temperature  of  the  hot  streak  was  twice  that 
of  the  surrounding  inlet  flow,  whereas  the  hot  streak  static  and 
stagnation  pressures  were  identical  to  the  free  stream.  The  hot 
streak  was  seeded  with  COi  and  the  path  of  the  hot  streak  de¬ 
termined  by  measuring  CO2  concentrations  at  various  locations 
within  the  turbine  stage  using  the  blade  surface  static  pressure 
taps.  In  the  numerical  simulation,  one  hot  streak  is  introduced 
to  the  inlet  of  every  stator  passage  in  the  form  of  a  hyperbolic 
tangent  (step-like)  temperature  profile.  A  hot  streak  tempera¬ 
ture  of  1.2  times  that  of  the  surrounding  inlet  flow  was  chosen 
for  this  investigation.  Shear  layer  instabilities,  caused  by  the 
step-like  temperature  profile,  developed  in  the  numerical  simu¬ 
lations  when  the  hot  streak  temperature  was  increased  beyond 
1.2  times  that  of  the  surrounding  flow.  The  predicted  tempera¬ 
ture  distributions  of  the  current  study  can  still  be  related  to  the 
experimental  data  using  the  temperature  coefficient  scaling  tech¬ 
nique  developed  by  Sharma  et  aJ  (1990).  In  the  1-stator/l-rotor 
simulation,  the  hot  streak  was  introduced  over  one  quarter  of  the 
stator  pitch  and  centered  at  midgap. 

For  the  1-stator/l-rotor  simulation,  the  stator  grid  system 
was  constructed  with  101  x  31  (streamwise  x  tangential)  grid 
points  in  the  O-grid  and  85  x  51  grid  points  in  the  H-grid.  The 
rotor  grid  system  was  constructed  with  101  x 31  grid  points  in  the 
O-grid  and  78x51  grid  points  in  the  H-grid.  A  total  of  14,575  grid 
points  were  used  in  this  simulation.  The  stator  airfoils  had  an 
average  y+  value  of  0.84,  while  the  rotor  airfoils  had  an  average 
y+  value  of  0.70.  Figure  1  illustrates  the  grid  topology  used  in 
the  numerical  simulations,  where  every  other  grid  point  in  the 
O-grid  has  been  omitted  for  clarity. 

The  numerical  simulations  performed  during  this  investiga¬ 
tion  were  computed  on  a  four  processor  Alliant  FX-80  mini¬ 
supercomputer.  Typical  calculations  required  .00191  seconds  per 
grid  point  per  iteration  computation  time.  In  the  current  study, 
approximately  six  cycles  at  2000  iterations  per  global  cycle  were 
needed  to  obtain  a  time-periodic  solution. 

In  the  experiment,  the  flow  coefficient  was  held  at  fa  = 
uc/U  =  .35.  The  subscript  ‘c’  refers  to  the  fact  that  the  flow 
coefficient  is  based  on  the  inlet  velocity  in  the  regions  outside 
the  hot  streak  (i.e.  cold  flow).  Within  the  hot  streak,  a  velocity 
increase  proportional  to  s/Th,/ T<x,  causes  a  significant  increase 


in  the  average  flow  coefficient,  0  =  u/U.  Thus,  in  the  numeri¬ 
cal  simulation,  where  a  hot  streak  is  included  in  every  passage, 
the  effect  of  this  velocity  increase  within  the  hot  streak  on  the 
flow  coefficient  must  be  taken  into  account.  In  addition,  previous 
experimental  data  (Dring  et  aL  1982)  and  numerical  simulations 
indicate  that  at  a  flow  coefficient  of  0  =  .35  the  flow  separates 
on  the  pressure  surface  of  the  rotor.  To  avoid  the  additional 
complications  that  flow  separation  may  present  in  interpreting 
the  numerical  results,  the  flow  coefficient  in  the  numerical  sim¬ 
ulation  was  increased  until  the  time-averaged  pressure  side  sep¬ 
aration  disappeared  (approximately  fa  —  .385).  The  difference 
in  the  flow  coefficients  equates  to  the  numerical  simulation  being 
performed  at  approximately  5  degrees  more  positive  incidence  at 
the  rotor  inlet  than  in  the  experiment.  This  increase  in  incidence 
corresponds  to  only  a  10%  shift  in  the  spanwise  location  (to¬ 
wards  the  tip)  at  which  the  numerical  simulation  was  performed, 
to  a  position  located  at  60%  span.  At  this  low  flow  coefficient,  a 
50  deg  variation  in  the  incidence  occurs  from  the  hub  to  the  tip 
of  the  rotor  airfoil.  This  can  be  shown  using  an  analytic  expres¬ 
sion  which  relates  the  rotor  inlet  relative  flow  angle  to  the  flow 
coefficient  (assuming  a  given  blade  geometry): 

cot  0  =  cot  a  -  <p~'  (6) 

where  a  is  the  absolute  flow  angle  at  the  stator  exit.  As  Figure  2 
illustrates,  the  rotor  inlet  flow  angle  changes  very  rapidly  with 
flow  coefficient  for  this  particular  geometry. 

In  order  to  understand  the  incidence  variation  at  the  rotor 
inlet  caused  by  the  hot  streak  and  stator  wakes,  a  study  was 
performed  to  determine  their  relevance.  By  taking  the  derivative 
of  eq.  (6),  an  equation  is  obtained  which  relates  the  change  in  the 
rotor  inlet  relative  flow  angle  to  the  change  in  the  flow  coefficient: 

A/3  =  A^/csc2#  (7) 

where  j}  is  the  mean  rotor  relative  inlet  flow  angle.  A0  is  the 
variation  in  the  flow  coefficient  due  to  a  hot  streak  or  stator 
wake,  and  A/3  is  the  change  in  the  rotor  inlet  relative  flow  angle. 


Figure  2:  Rotor  inlet  relative  flow  angle  vs  Sow  coefficient 
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Figure  3:  Rotor  inlet  flow  angle  and  static  temperature  distribu-  Figure  4:  Predicted  pressure  coefficient  distribution  for  the  stator 

tions 


Equation  (7)  is  useful  in  estimating  the  change  in  the  rotor  inlet 
relative  flow  angle  as  the  rotor  passes  through  the  hot  streak.  For 
the  present  numerical  simulation,  eq.  (7)  predicts  a  4.12  deg  vari¬ 
ation  in  the  rotor  inlet  relative  flow  angle  due  to  the  hot  streak. 
Figure  3  shows  the  numerically  predicted  rotor  inlet  relative  flow 
angle  and  static  temperature  for  one  cycle,  where  a  cycle  cor¬ 
responds  to  the  rotor  blade  rotating  through  an  angle  of  2x/N 
where  N  is  the  number  of  stator  blades  (i.e.  N=28  for  the  present 
simulations).  The  stator  wake  has  little  effect  on  the  rotor  inlet 
static  temperature,  but  has  a  significant  impact  on  the  relative 
flow  angle.  The  wake  causes  the  rotor  inlet  relative  flow  angle  to 
vary  by  approximately  4.9  deg  in  the  direction  of  more  negative 
incidence  and  affects  the  flow  angle  for  nearly  50%  of  the  cycle. 
The  hot  streak  is  seen  to  cause  a  variation  in  both  the  static  tem¬ 
perature  and  relative  inlet  flow  angle.  The  hot  streak  causes  the 
rotor  inlet  relative  flow  angle  to  vary  by  approximately  3.9  deg 
(compared  to  4.12  deg  calculated  above)  in  the  direction  of  more 
positive  incidence  and  affects  the  flow  angle  for  about  50%  of  the 
cycle.  Thus,  through  the  course  of  one  cycle  the  rotor  inlet  rel¬ 
ative  flow  angle  varies  by  almost  10 deg,  although  the  incidence 
variation  in  the  rotor  caused  by  the  stator  wake  is  nearly  can¬ 
celled  by  that  corresponding  to  the  hot  streak  for  time-averaged 
flow.  Figure  3  also  shows  that  the  flow  angle  variation  due  to  the 
hot  streak  lags  the  temperature  increase  due  to  the  hot  streak. 


Figure  4  illustrates  the  predicted  stator  surface  time-averaged 
pressure  coefficient  distribution.  Unfortunately,  no  experimental 
time-averaged  pressure  data  was  available  for  the  stator  airfoils. 
The  time-averaged  pressure  coefficient  is  defined  as: 


C  —  ~  Fn ) 


PxlP 


(8) 


where  P«*»  is  the  local  time-averaged  pressure,  P,x  is  the  inlet 
total  pressure,  p\  is  the  average  inlet  free  stream  density,  and  V 
is  the  rotor  velocity.  Figure  5  shows  the  predicted  and  experi¬ 
mental  (Dring  et  a i,  1982)  rotor  surface  time-averaged  pressure 
coefficient  distributions.  While  the  experimental  data  of  Dring  et 
of  ( 1982)  is  for  rotor/stator  flow  without  a  hot  streak,  past  sim¬ 
ulations  (Doroey  ef  of,  1990,  Rai  and  Dring,  1990)  have  shown 
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Figure  5:  Predicted  and  experimental  pressure  coefficient  distri¬ 
butions  for  the  rotor 

that  for  the  present  type  of  simulation,  the  pressure  and  temper¬ 
ature  fields  are  decoupled.  Excellent  agreement  exists  between 
the  predicted  results  and  the  experimental  data  for  the  rotor, 
except  near  the  suction  surface  leading  edge  where  the  numeri¬ 
cal  simulation  predicts  a  strong  overspeed.  This  discrepancy  is 
probably  due  to  the  slight  difference  (1  deg)  in  the  time- aver  aged 
incidence  to  the  rotor  between  the  numerical  simulation  and  the 
experiment  (Dring  et  al  1982). 

The  migration  path  of  the  numerical  hot  streak  is  displayed  in 
Fig.  6,  which  shows  static  temperature  contours  at  one  instant 
in  time  during  the  course  of  a  cycle.  The  hot  streak  is  undis¬ 
turbed  as  it  migrates  through  the  stator  passage,  except  that  the 
width  of  the  hot  streak  decreases  due  to  flow  acceleration.  The 
geometry  of  the  rotor  blade  is  such  that  the  hot  streak  first  im¬ 
pinges  upon  the  pressure  surface,  then  wraps  around  the  leading 
edge  and  moves  along  the  suction  surface.  The  hot  fluid  on  the 
suction  surface  is  quickly  convected  downstream,  while  the  hot 
fluid  near  the  pressure  surface  remains  in  the  rotor  passage  for 
some  time  after  the  rotor  passes  through  the  hot  streak  The  hot 
streak  assumes  a  V  shape  (due  to  the  difference  in  convection 
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Figure  6:  Static  temperature  contours  for  the  LSRR  turbine  stage 


speeds  between  the  pressure  and  suctions  sides  of  the  passage) 
as  it  continues  to  move  downstream  through  the  rotor  passage, 
before  breaking  up  and  leaving  the  turbine  in  the  form  of  discrete 
temperature  eddies. 

Figure  7  illustrates  the  experimental  time-averaged  C02  con¬ 
centration  contours  for  the  rotor  surface.  The  movement  of  the 
hot  streak  fluid  from  leading  edge  to  the  trailing  edge  of  the 
pressure  surface  at  approximately  constant  spanwise  locations 
indicates  the  two-dimensional  nature  of  the  hot  streak.  Con¬ 
spicuously  absent  is  the  radial  migration  of  the  hot  streak  on 
the  pressure  surface  observed  both  experimentally  (Butler  et  al 
1989)  and  in  previous  numerical  simulations  (Dornev  et  al.  1990. 
Takahashi  and  Ni.  1990)  for  a  circular  hot  streak.  Figure  8  com¬ 
pares  the  predicted  time-averaged  temperature  coefficient  distri¬ 
bution  along  the  surface  of  the  rotor  to  the  experimental  data 
taken  at  midspan.  The  temperature  coefficient.  Cj ■  is  defined  as 
(Takahashi  and  Ni.  1990,  Sharma  et  al.  1990): 


Ct  = 


T-T, 


(9) 


where  T  is  the  local  time-averaged  temperature  and  T^grle  is 
the  time-averaged  temperature  at  the  rotor  leading  edge.  Re¬ 
sults  of  the  numerical  simulation  show  fair  agreement  with  the 


Figure  8:  Predicted  and  experimental  rotor  surface  temperature 

-  Ct 

planar  hot  streak  experimental  data.  The  numerical  -imulat;  u, 
underpredicts  the  temperature  near  the  pressure  -urt'ace  bail¬ 
ing  edge,  and  also  predicts  a  rapid  temperature  fail  on  un  ’tie 
suction  surface  which  is  not  observed  in  the  experimental  .lata 
Figure  9  shows  the  time-averaged  temperature  contour-  fur 

entire  rotor  passage.  The  accumulation  of  hot  streak  fluid  on  the 
pressure  side  of  the  passage  is  evident.  A  close-up  of  the  rui.-r 
pressure  surface  leading  edge  region  (see  Figure  10)  shows  that  a 
region  of  cooler  fluid  is  trapped  between  the  hot  streak  fluid  and 
the  rotor  surface.  Similar  results  were  reported  in  Dorney  et  al 
(1990).  for  a  two-dimensional  hot  streak  simulation  at  different 
operating  conditions.  Figure  11.  which  shows  time-averaged  tern 
perature  profiles  at  three  locations  along  the  pressure  surface  it 
the  rotor,  also  indicates  the  presence  of  a  cooler  region  near  'he 
surface  at  the  leading  edge.  Time-averaged  boundary  laser 


Figure  7:  Experimental  time-averaged  C02  contours  for  the 
LSRR  rotor 


Figure  9:  Time-averaged  temperature  contours  for  th-  - 
sage 


Figure  10:  Time-averaged  temperature  contours  near  the  rotor 
pressure  surface  leading  edge 

locations  are  also  included  in  Figure  11.  At  255?  axial  chord  the 
maximum  temperature  in  the  boundary  layer  occurs  off  the  sur¬ 
face.  approximately  one-third  the  distance  to  the  boundary  layer 
edge.  Similarly,  at  50%  axial  chord  the  maximum  temperature 
in  the  boundary  layer  occurs  well  above  the  rotor  surface.  At 
759?  axial  chord,  however,  the  time-averaged  boundary  layer  is 
noticeably  thinner  and  the  maximum  time-averaged  temperature 
occurs  at  the  rotor  surface.  A  general  observation  drawn  from 
this  investigation  is  that  as  the  boundary  layer  becomes  thin¬ 
ner.  the  maximum  temperature  in  a  given  boundary  layer  profile 
moves  closer  to  the  airfoil  surface. 

The  differences  between  the  experimental  data  and  the  pre¬ 
dicted  surface  temperatures  shown  in  Fig.  8  may  be  due  to  one 
of  several  factors.  As  shown  in  Dorney  et  al  (1990),  Takahashi 
and  N'i  ( 1990),  and  Sharma  et  al  ( 1990)  a  significant  effect  on  the 
surface  temperature  exists  (for  certain  hot  streak  shapes)  due  to 
the  three-dimensionality  of  the  rotor  surface  boundary  layers  and 
the  influence  of  secondary  flows.  A  second  factor  which  may  also 
contribute  to  the  discrepancies  between  the  experimental  data 


Figure  11:  Predicted  time-averaged  temperature  profiles 


and  the  two-dimensional  simulation  is  the  assumed  correlation 
between  the  predicted  surface  temperature  and  the  experimen¬ 
tally  measured  C0 1  concentrations.  In  the  experiment,  the  CO : 
concentrations  along  the  surface  of  the  rotor  were  determined  by 
drawing  samples  of  the  gas  in  through  static  pressure  taps.  Tr.e 
suction  force  may  cause  CO2  gas  from  well  above  the  airfoil  sur¬ 
face  to  be  included  in  the  sample.  To  test  this  hvpothesis.  'lie 
time-averaged  temperature  coefficient  was  redefined  as: 


C  t 


tu  -  r. 


IU  ■ 


where  Ti,;  is  the  time-averaged  temperature,  area-averaged  oier 
a  given  boundary  layer  profile.  Figure  12  compares  the  predicted 
results  using  the  modified  temperature  coefficient  definition  with 
the  experimental  data.  Excellent  agreement  now  exists  between 
the  predicted  and  experimental  results.  A  comparison  of  F.g- 
ures  8  and  12  reveals  that  the  predicted  surface  tempera- ur*- 
is  approximately  equal  to  the  area-averaged  temperature  in  the 
boundary  layer,  except  in  the  leading  edge  region.  Figure  '.2 
shows  that  two-dimensional  simulations  can  be  used  to  provide 
guidance  as  to  when  and  where  rotor  surface  temperatures  ex 
ceed  allowable  Emits,  provided  an  averaged  temperature  aero-? 
the  boundary  layer  is  examined.  This  is  an  important  finding 
which  shows  that  costly  and  time  consuming  three-dimensional 
simulations  of  hot  streak  migration  may  not  need  to  be  executed 
to  determine  if  rotor  pressure  surface  burning  is  likeiy .  Further 
investigation  is  required  to  determine  if  this  manner  of  examining 
average  viscous  layer  temperatures  is  reliable  and  applicable 
other  flow  conditions  and  geometries. 


Figure  12:  Predicted  and  experimental  rotor  su:n> 
t  ures  -  Ct 


CONCLUSIONS 

A  numerical  investigation  has  been  performed  to  determine  if 
two-dimensional  stator/rotor/hot  streak  interaction  simulations 
can  provide  reliable  warnings  to  turbine  designers  as  to  when  and 
where  rotor  surface  temperatures  exceed  allowable  limits.  Dur¬ 
ing  this  investigation,  the  results  from  an  unsteady  1-stator '1- 
rotor/l-hot  streak  Navier-Stokes  simulation  were  compared  with 
new  experimental  data  taken  from  a  turbine  stage  v.-hich  included 
a  planar  hot  streak  introduced  at  the  stator  inlet.  The  predicted 
and  experimental  rotor  surface  temperature  distributions  show- 
close  agreement,  provided  the  proper  correlation  is  made  be- 
ttveen  the  experimental  CO 2  concentrations  and  the  numerically 
predicted  temperatures.  Thus,  for  certain  hot  streak  shapes,  two- 
dimensional  numerical  procedures  can  be  used  as  accurate  and 
efficient  tools  by  designers  attempting  to  optimize  cooling  sys¬ 
tems. 

Future  work  will  focus  on  developing  improved  techniques  for 
comparing  numerically  predicted  temperature  distributions  with 
experimental  data.  Based  on  the  agreement  between  the  current 
numerical  results  with  the  available  experimental  data,  another 
investigation  will  be  performed  which  uses  the  two-dimensional 
Navier-Stokes  procedure  to  identify  the  effects  of  rotor  inlet  rel¬ 
ative  flow  angle  and  boundary  layer  thickness  on  hot  streak  mi¬ 
gration  and  rotor  pressure  surface  heating. 
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Summary 

A  three-dimensional  Navier-Stokes  analysis  of  heat  transfer  and 
aerodynamic  performance  is  presented  for  a  low  speed  linear 
turbine  cascade.  The  numerical  approach  used  in  this  anal¬ 
ysis  consists  of  an  alternate-direction,  implicit,  approximate- 
factorization,  time-marching  technique.  An  objective  of  this  in¬ 
vestigation  has  been  to  establish  the  computational  grid  density 
requirements  necessary  to  accurately  predict  blade  surface  and 
endwall  heat  transfer,  as  well  as  the  exit  plane  aerodynamic  to¬ 
tal  pressure  loss  and  flow  angle  distributions.  In  addition,  a  study 
has  been  performed  to  determine  a  viable  implementation  strat¬ 
egy  for  the  three-dimensional  modelling  of  transition  and  turbu¬ 
lence  in  the  turbine  blade  passage.  Results  are  presented  which 
demonstrate  that  the  present  procedure  can  accurately  predict 
3-D  turbine  blade  heat  transfer  as  well  as  the  absolute  level  and 
spanwise  distribution  of  aerodynamic  performance  quantities. 

List  of  Symbols 
a  -  Speed  of  sound 

c,  -  Specific  heat  ratio  at  constant  preaure 
e  -  Specific  energy 

e,  -  Total  energy 

k  -  Thermal  conductivity 

L  -  Reference  length  (2.54  cm) 

M  -  Mach  number 

P  -  Static  pressure 

Pt  -  Prandtl  number 

Pj-  ■  Stagnation  pressure 

f  -  Heat  flux 

He  -  Free  stream  inlet  reference  Reynolds  number 
T  -  Static  temperature 

u  -  x  component  of  velocity 

•  -  y  component  of  velocity 

w  •  z  component  of  velocity 

X  -  Second  coefficient  of  viscosity 

M  -  First  coefficient  of  viseoaty 

f  -  Density 

r  -  Shear  stress 

Subscripts 

t  -  Inviscid 

L  -  Laminar  quantity 

t  -  Stagnation  quantity 

T  -  Turbulent  quantity 

v  -  Viscous 

s,f,x  •  First  derivative  with  respect  to  x,  y  or  z 

sz,  yy,  xx  -  Second  derivative  with  respect  to  x,  y  or  t 

1  -  Inlet  quantity 

2  -  Exit  quantity 
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Introduction 

The  accurate  prediction  of  three-dimensional  turbomachinery 
blade  row  heat  transfer  and  aerodynamic  performance,  including 
total  pressure  loss  and  flow  turning,  remains  a  challenging  task 
for  most  Navier-Stokes  numerical  procedures.  The  very  large 
computational  grid  densities  necessary  to  accurately  resolve  the 
viscous  flow  physics  and  properly  convect  the  vorticity  require 
computational  resources  which  far  exceed  those  available  to  most 
engineering  facilities.  Thus,  in  most  instances,  it  has  not  been 
possible  to  produce  grid-independent  three-dimensional  flow  so¬ 
lutions  for  turbomachinery  blade  rows  which  ensure  the  accuracy 
required  for  reliable  predictions. 

The  proper  modelling  of  transition  and  turbulence  within  tur¬ 
bomachinery  blade  passages  is  also  often  compromised  in  order 
to  reduce  computational  times  to  acceptable  levels.  The  effects  of 
turbulence  in  the  flow  are  quite  often  modelled  with  an  algebraic 
eddy  viscosity  correlation  which  has  its  basis  in  two-dimensional 
boundary  layer  data,  even  though  the  flow  in  a  blade  row  passage 
has  strong  secondary  and  occasional  reversed  flows.  The  mod¬ 
elling  of  transition  in  most  Navier-Stokes  approaches  is  practi¬ 
cally  non-existent  due  to  the  lack  of  reliable  correlations  for  three- 
dimensional  internal  blade  row  flows  and  due  to  the  numerical 
instabilities  caused  by  the  unsteady  movement  of  the  transition 
region  as  the  flow  field  develops.  Unfortunately,  it  is  unlikely 
that  improvements  in  turbulence  and  transition  modelling  will 
improve  beyond  the  current  state-of-the-art  for  three-dimensional 
turbomachinery  flows  until  Navier-Stokes  procedures  can  afford¬ 
ably  and  reliably  produce  grid-independent  solutions.  Only  then 
will  researchers  be  able  to  separate  out  the  turbulence  modelling 
effects  on  the  flow  from  those  caused  by  truncation  error  and 
numerical  viscosity.  * 

Several  investigations  have  been  made  recently  in  which  three- 
dimensional  Navier-Stokes  procedures  have  been  used  to  pre¬ 
dict  the  heat  transfer  and/or  aerodynamic  performance  of  the 
Langston  linear  turbine  cascade  [1,  2J.  Among  these  numerical 
investigations.  Chan  and  Sheedy  [3]  have  used  a  pressure  cor¬ 
rection  method  with  a  two  equation  (t  -  t)  turbulence  model 
to  predict  the  losses  in  the  Langston  cascade  for  a  thick  (3.3 
cm)  inlet  boundary  layer.  Their  predicted  results  indicate  the 
numerical  procedure  is  very  sensitive  to  spanwise  computational 
grid  density.  Even  the  use  of  comparatively  fine  computational 
grids,  however,  was  not  adequate  to  achieve  the  flow  field  res¬ 
olution  needed  to  accurately  predict  aerodynamic  preformance 
quantities.  Hah  [4,  5]  has  also  used  a  pressure  correction  tech¬ 
nique  with  a  two  equation  ( t  -  <)  turbulence  model  to  predict 
endwall  Stanton  numbers  [4]  and  losses  as  a  function  of  axial  dis¬ 
tance  [5)  for  the  Langston  cascade  with  a  thick  inlet  bonndary 
layer.  The  predicted  endwall  Stanton  numbers  show  fair  agree¬ 
ment  with  the  experimental  data,  but  it  is  unclear  if  a  sufficiently 
fine  computational  grid  was  used  to  attain  grid  independent  solu¬ 
tions.  Moore  [6,  7]  has  used  a  pressure  correction  method  with  an 
algebraic  turbulence  model  to  predict  losses  and  secondary  Sows 
for  the  Langston  cascade  with  a  thick  inlet  boundary  layer.  Lee 
and  Knight  [8]  solved  the  thin  layer  Navier-Stokes  equations  us¬ 
ing  an  approximate  factorization  technique  and  the  two  equation 
( q  -  w)  turbulence  model  to  study  heat  transfer  for  the  Langston 
cascade  with  a  thick  inlet  boundary  layer.  A  portion  of  their 
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investigation  focused  on  the  effects  of  different  grid  topologies 
(H-  and  0-H  grids)  on  the  predicted  heat  transfer  quantities. 
In  their  study,  it  was  determined  that  grid  topology,  as  well  as 
computational  grid  density,  can  impact  the  predicted  solutions. 

The  objective  of  the  current  investigation  has  been  to  pro¬ 
vide  guidance  as  to  the  computational  grid  density  requirements 
which  are  necessary  to  produce  reliable  and  accurate  Navier- 
Stokes  predictions  of  turbomachinery  blade  row  flows.  Simu¬ 
lations  of  transitional  flow  through  the  Langston  cascade  (l,  2] 
with  a  thin  inlet  boundary  layer  (0.41  cm)  have  been  performed 
where  the  focus  has  been  to  conduct  a  computational  grid  refine¬ 
ment  study  to  determine  the  computational  grid  density  neces¬ 
sary  to  properly  resolve  surface  heat  transfer  and  aerodynamic 
performance  quantities.  As  part  of  this  investigation,  a  study 
was  performed  to  determine  a  viable  strategy  for  implement¬ 
ing  an  algebraic  turbulence  model  for  three-dimensional  inter¬ 
nal  blade  row  flows  and  to  determine  the  limitations  of  using 
an  algebraic  turbulence  model  on  the  accuracy  of  the  exit  plane 
distributions  of  total  pressure  loss  and  flow  angle.  Also,  a  brief 
study  was  conducted  to  demonstrate  the  importance  of  modelling 
laminar/turbuient  transition  for  these  flows. 


Numerical  Integration  Procedure 


The  governing  equations  considered  in  this  study  are  the  time 
dependent,  three-dimensional  Navier-Stokes  equations: 


V*  +  (f.  +  K)z  +  (C,  -I-  Gu)y  +  (ff.  +  H,)t  =0  (1) 
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For  the  present  application,  the  second  coefficient  of  viscosity 
is  calculated  using  Stokes'  hypothesis,  A  =  -2/3/1.  The  equations 
of  motion  are  completed  by  the  perfect  gas  law. 

The  viscous  fluxes  are  simplified  by  incorporating  the  thin 
layer  assumption  [9].  In  the  current  study,  viscous  terms  are 
retained  in  the  direction  normal  to  the  hub  surface  (z-direction) 
and  in  the  direction  normal  to  the  blade  surface  (y-direction). 
To  extend  the  equations  of  motion  to  turbulent  flows,  an  eddy 
viscosity  formulation  is  used.  Thus,  the  effective  viscosity  and 
effective  thermal  conductivity  can  be  defined  as: 


P  =  PL  +  PT 


(7) 


The  turbulent  viscosity,  pj,  is  calculated  using  the  Baldwin- 
Lomax  [9]  algebraic  turbulence  model.  The  specific  techniques 
used  to  implement  this  model  for  three-dimensional  flows  in  tur- 
bomachinery  blade  rows  will  be  described  below. 

The  numerical  procedure  for  the  three-dimensional  analysis 
consists  of  a  time  marching,  implicit,  third-order  spatially  accu¬ 
rate.  upwind,  finite  difference  scheme.  The  inviscid  fluxes  are 
discretized  using  a  combination  of  central,  forward,  and  back¬ 
ward  differences  depending  on  the  local  eigenvalues  of  the  flux 
Jacobians  according  to  the  scheme  developed  by  Roe  [lOj.  The 
viscous  fluxes  are  calculated  using  standard  central  differences. 
The  alternate  direction,  approximate-factorization  technique  of 
Beam  and  Warming  [11]  is  used  to  compute  the  time  rate  changes 
in  the  primary  variables.  An  inner  Newton  iteration  can  be  used 
to  increase  stability  and  reduce  linearization  errors.  For  all  cases 
investigated  in  this  study,  one  Newton  iteration  was  performed  at 
each  time  step.  Further  details  of  the  three-dimensional  numeri¬ 
cal  technique  can  be  found  in  Refs.  [12,  13]. 


Grid  Generation  and  Geometry 


The  three-dimensional  Navier-Stokes  analysis  uses  zonal  grids  to 
discretize  the  turbine  flow  field.  A  combination  of  0-  and  H-grid 
sections  are  generated  at  constant  radial  spanwise  locations  in  the 
blade-to-blade  direction  extending  upstream  of  the  airfoil  lead¬ 
ing  edge  to  downstream  of  the  airfoil  trailing  edge.  Algebraically 
generated  H-grids  are  used  in  the  regions  upstream  of  the  lead¬ 
ing  edge,  downstream  of  the  trailing  edge  and  in  the  inter-blade 
region.  The  O-grid.  which  is  body-fitted  to  the  surface  of  the  air¬ 
foil  and  generated  using  an  elliptic  equation  solution  procedure,  is 
used  to  properly  resolve  the  viscous  flow  in  the  blade  passages  and 
to  easily  apply  the  algebraic  turbulence  model.  Computational 
grid  lines  within  the  O-grid  are  stretched  in  the  blade-normal  di¬ 
rection  with  a  fine  grid  spacing  at  the  wall.  The  combined  H-  and 
0-  overlaid  grid  sections  are  stretched  in  the  spanwise  direction 
away  from  the  hub  and  shroud  regions  with  a  fine  grid  spacing 
located  adjacent  to  the  hub  and  shroud. 


i 
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Boundary  Conditions 

The  theory  of  characteristics  is  used  to  determine  the  boundary 
conditions  at  the  passage  inlet  and  exit.  For  subsonic  inlet  flow, 
the  total  pressure,  o  and  tu  velocity  components  and  the  down¬ 
stream  running  Riemann  invariant,  R i  =  u  +  are  specified 
while  the  upstream  running  Riemann  invariant.  Rj  =  u  -  is 
extrapolated  from  the  interior  of  the  computational  domain.  For 
subsonic  outflow  the  pressure  ratio.  Pi/Pt\,  is  specified  while  the 
v  and  ui  velocity  components,  entropy,  and  the  downstream  run¬ 
ning  Riemann  invariant  are  extrapolated  from  the  interior  of  the 
computational  domain.  Periodicity  is  enforced  along  the  outer 
boundaries  of  the  H-grids  in  the  circumferential  (8)  direction. 


the  parallel  distance  is  measured  from  the  blade  surface  to  the 
grid  point  and  the  normal  distance  is  measured  from  the  endwall 
to  the  grid  point.  The  eddy  viscosity  is  then  based  on  the  flow 
variables  along  a  computational  grid  line  from  the  endwall  to  the 
grid  point.  Two  significant  problems  arise  from  this  particular 
three-dimensional  implementation  of  the  Baidwin-Lomax  turbu¬ 
lence  model.  First,  the  eddy  viscosity  distribution  in  the  cor¬ 
ner  regions  is  discontinuous  across  the  J  =  K  computational  lines 
which  can  cause  large  gradients  to  occur  in  the  velocity  field. 
Secondly,  this  particular  blending  is  highly  dependent  upon  the 
computational  grid  density  and  stretching  in  both  computational 
directions.  It  was  determined,  based  upon  numerical  simulations, 
that  the  steady-state  solution  in  the  blade/endwall  region  was  ex¬ 
tremely  sensitive  to  changes  in  the  computational  grid  structure. 


Dirichlet  conditions,  in  which  the  time  rate  change  in  the 
vector  V  of  Eq.  (2)  is  set  to  zero,  are  imposed  at  the  overlaid 
boundaries  of  the  0-  and  H-grids.  The  flow  variables  of  U  at 
zonal  boundaries  are  explicity  updated  after  each  time  step  by 
interpolating  values  from  the  adjacent  grid.  Because  of  the  ex¬ 
plicit  application  of  the  zonal  boundary  conditions,  large  time 
steps  may  necessitate  the  use  of  more  than  one  Newton  itera¬ 
tion.  The  zonal  boundary  conditions  are  non-conservative,  but 
for  subsonic  flow  this  should  not  affect  the  accuracy  of  the  final 
flow-  solution. 

No-slip  boundary  conditions  are  enforced  at  the  hub  and 
shroud  endwalls  of  the  turbine  passage  and  along  the  surface  of 
the  airfoil.  In  addition,  a  specified  heat  flux  or  wall  temperature 
distribution  is  held  constant  in  time  along  the  solid  surfaces.  In 
the  present  investigation  a  specified  heat  flux  given  by 

.OT 

Vwall  —  *  toy 

is  prescribed  over  all  solid  surfaces  at  a  level  consistent  with  the 
experimental  conditions.  Upon  descretizing  Eq.  8  using  a  first  or¬ 
der  backward  difference,  a  Neumann  boundary  coaditioa  arises 
which  is  held  during  the  implicit  sweep  of  lines  normal  to  the 
surface.  The  implicit  update  of  the  solid  boundary  flow  variables 
is  followed  by  a  second  order  accurate  explicit  update.  Further 
information  describing  the  implementation  of  the  boundary  con¬ 
ditions  can  be  found  in  Refs.  [12,  13,  14], 

Turbulence/Transition  Model 

The  Bald  win- Lomax  [9]  two-layer  algebraic  turbulence  model  is 
used  to  represent  the  turbulent  eddy  viscosity  in  the  flow  field. 
This  turbulence  model  is  based  upon  two-dimensional  boundary 
layer  data,  and  as  such,  is  not  well  suited  for  comer  flows  such 
as  those  at  the  blade/endwall  juncture.  Previously,  the  treat¬ 
ment  used  to  implement  this  turbulence  model  in  the  comer  re¬ 
gions  was  the  technique  proposed  by  Hung  and  Buning  [15,  16], 
In  this  technique,  the  Baidwin-Lomax  turbulence  model  is  com¬ 
puted  separately  for  each  endwall  and  the  blade  surface.  The 
mixing  length  in  the  comer  region  depends  on  the  computational 
indicies  of  a  given  node.  For  instance,  consider  the  cue  when  the 
inconstant  computational  lines  run  normal  to  the  blade  surface 
and  the  K=constant  Unes  run  normal  to  the  endwall.  For  any 
computational  node  whose  J-wise  index  is  leu  than  its  K-wise 
index,  the  ‘normal’  distance  is  defined  u  the  distance  from  the 
blade  surface  to  the  grid  point  and  the  ‘parallel’  distance  is  de¬ 
fined  as  the  distance  from  the  endwall  to  the  grid  point.  The 
mixing  length  for  the  inner  region  of  the  boundary  layer  is  then 
calculated  as 

I  =  2sn/(s  +  n  +  ^/js^-fn*))  (9) 

where  s  is  the  parallel  distance  and  n  is  the  normal  distance.  The 
eddy  viscosity  is  then  based  on  the  flow  variables  along  a  com¬ 
putational  grid  line  from  the  airfoil  surface  (normal  direction)  to 
the  grid  point  under  consideration.  Likewiu,  for  any  computa¬ 
tional  node  whose  J-wise  index  is  greater  than  its  K-wise  index. 


In  the  present  investigation,  a  blending  function  has  been 
used  to  smoothly  vary  the  eddy  viscosity  distribution  between 
the  blade  and  endwall.  Separate  eddy  viscosity  distributions  are 
computed  for  the  blade  and  endwall  surfaces  along  the  computa¬ 
tional  lines  which  run  normal  to  each  surface  respectively.  The 
eddy  viscosity  in  the  corner  flow  regions  is  then  computed  based 
upon  the  following  blending  function  according  to  the  work  of 
Vatsa  and  Wedan  [17]: 


„  _  jLiHi  t  ^**r- 
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(10) 


where  dt,  is  the  distance  from  the  blade  surface  to  a  given  node, 
dew  is  the  distance  from  the  endwall  surface  to  the  node,  and 
pr4,  P7,„  are  the  eddy  viscosities  computed  from  the  blade  and 
endwall  flows,  respectively.  This  type  of  blending  greatly  reduces 
the  dependence  of  the  eddy  viscosity  distribution  on  the  compu¬ 
tational  grid  and  creates  a  smooth  eddy  viscosity  distribution  in 
the  comer  regions. 

As  will  be  shown  below,  results  from  a  number  of  computa¬ 
tions  demonstrate  the  importance  of  correctly  modelling  transi¬ 
tion  of  the  flow  from  laminar  to  turbulent  on  the  blade  surface. 
In  the  present  investigation  in  which  turbine  blade  heat  trans¬ 
fer  is  being  predicted,  the  modelling  of  transition  is  essential  for 
the  accurate  prediction  of  Stanton  number  distributions.  In  this 
study,  transition  was  modelled  using  a  prescribed  location  for  the 
onset  of  transition.  The  onset  of  transition  was  varied  as  a  func¬ 
tion  of  span,  with  the  endwalls  assumed  to  be  fully  turbulent 
and  mid-span  transition  beginning  at  the  20%  axial  chord  loca¬ 
tion  on  the  suction  surface  according  to  the  experimental  data. 
Transition  was  assumed  to  occur  as  an  instantaneous  process.  As 
will  be  discussed  below,  a  means  of  accurately  predicting  transi¬ 
tion  especially  in  turbine  blade  rows  clearly  must  be  developed 
in  order  to  reliably  and  accurately  predict  heat  transfer. 


Results 


A  set  of  viscous  three-dimensional  numerical  simulations  of  flow 
through  the  Langston  low  speed  turbine  cascade  with  a  thin  inlet 
boundary  layer  (0.41  cm),  including  heat  transfer  and  transition, 
have  been  performed  and  the  predicted  results  compared  to  ex¬ 
perimental  data  [1.  2],  The  heat  flux  values  at  the  endwall  and 
airfoil  surfaces  were  specified  to  be  the  same  as  in  Ref.  [2].  The  in¬ 
let  Mach  number  was  set  to  Afj  =  0.0917  and  tbe  inlet  flow  angle 
was  specified  to  be  0\  =  44.7  degrees.  The  inlet  Reynolds  num¬ 
ber  was  set  to  Re:  =  25590  per  cm  and  the  ratio  of  the  exit  static 
pressure  to  the  inlet  total  pressure  was  fixed  at  Fj//*ri  *  0.9847 
based  on  the  experimental  data. 

A  computational  grid  refinement  study  has  been  performed 
to  establish  the  grid  density  required  to  produce  essentially  grid- 
independent  solutions  in  both  the  blsde- to- blade  and  spanwise 
directions.  A  combination  of  two-  and  three-dimensional  simula¬ 
tions  were  performed  to  establish  the  reqoired  grid  density  and 
wall  spacing  in  the  blade-to-blade  direction  to  accurately  predict 
turbine  blade  heat  transfer.  As  a  result  of  the  two-dimensional 
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Case 

Grid  Density 
stream  xtangx  span 

•^Vbladc 

1 

Aj/cnJuaii 

'IT » 

Profile 

APT  PT ; 

Total 

APT/PT: 

Profile 

3: 

Total 

32 

Trns 

i 

H-70  x  31  x  25 
0-101  x  31  x  25 

0.00100 

(4.260) 

0.00300 

1 7  430i 

0.20*46 

0.19065 

27.96 

27  34 

yes 

2 

H-70  x  31  x  25 
0-101  x  31  x  25 

0.00010 

1 0.3*4 ) 

0.00100 

.  2.920  i 

0.16102 

0.1*373 

27.49 

27.36 

no 

3 

H-70  x  31  x  25 
0-101  x  31  x  25 

0.00010 

10.379) 

0.00100 

1 2.761  p 

0.16100 

0 1*372 

27.49 

27  5* 

yp> 

4 

H-70  x  31  x  51 
0-101  x  31  x  51 

0.00010 

10.360) 

0.00010 

<0.266) 

0.14575 

0.17271 

27.67 

27.36 

yes 

5 

H-70  x  31  x  71 
0-101  x  31  x  71 

0.00010 

10.351) 

0.00010 

1 0.266) 

0  13696 

0.17075 

27.73 

27  26 

>ps 

6 

H-70  x  31  x  91 
0-101  x  31  x  91 

0.00010 

0.00010 

0.13051 

0.17064 

27.75 

27.25 

yes 

"D 

0-131  x  61 

H-71  x  21 

0.00002 

10.0591 

— 

0.12030 

26.56 

yes 

r..\PT 

— 

- 

0.13000 

0.17437 

26.60 

26.32 

— 

Table  1:  Grid  refinement  study 


simulations,  it  was  determined  that  a  computational  grid  consist¬ 
ing  of  70  x  21  grid  points  in  the  H-grid  and  101  x  61  grid  points 
in  the  O-grid  with  an  average  wall  spacing  of  .00002  (if*  =  .059) 
produced  grid-independent  solutions.  Once  this  was  established, 
a  similar  grid  topology  and  density  were  used  in  the  blade-to- 
blade  direction  for  the  three-dimensional  simulations.  The  wall 
spacing  adjacent  to  the  airfoil  surface  was  then  refined  during  a 
series  of  three-dimensional  simulations  to  insure  the  accuracy  of 
the  three-dimensional  blade-to-blade  solution.  Finally,  the  span- 
wise  computational  grid  was  refined  both  in  terms  of  endwall 
spacing  and  density  until  the  exit  total  pressure  loss  and  exit 
flow  angle  distributions  remained  unchanged  and  exhibited  close 
agreement  with  the  experimental  data. 

A  summary  of  the  cases  executed  during  this  investigation  are 
given  in  Table  1.  The  columns  labelled  'Total'  in  Table  1  indicate 
area-averaged  values  in  the  spanwise  direction,  while  the  columns 
labelled  Profile'  refer  to  mid-span  values.  The  area-averaged  val¬ 
ues  were  calculated  using  trapezoidal  integration  and  only  the 
computational  grid  points  at  the  experimental  data  locations 
were  included  in  the  integration  (a  total  of  13  locations  for  the 
half-span)  in  order  to  be  consistent  with  the  experimental  area- 
averaged  exit  loss  and  flow  angle  values.  The  column  labelled 
Tims'  indicates  whether  or  not  transition  waj  modelled.  Since 
the  geometry  was  symmetric  about  the  mid-span,  the  numerical 
simulations  were  performed  for  only  the  half-span  and  a  symme¬ 
try  boundary  condition  was  used  along  the  symmetry  plane.  The 
spanwise  grid  densities  given  in  Table  1  are  for  the  half-span  ge¬ 
ometry.  Figure  1.  which  illustrates  the  Total'  exit  flow  angle  and 
total  pressure  toss  (tee  Eq.  13)  values  as  a  function  of  grid  point 
density,  demonstrates  that  a  half-span  grid-independent  solution 
was  obtained  with  approximately  450.000  grid  points.  Except 
where  noted,  the  predicted  results  shown  below  are  for  this  fine 
grid  (Case  6)  solution- 


Figure  1:  Exit  flow  angle  and  total  pressure  loss  as  a  function  of 
grid  point  density 


Figure  2:  Blade-to-blade  computational  grid  for  Langston  cas¬ 
cade 

Figure  2  illustrates  an  example  of  the  blade-to-biade  com¬ 
putational  grid  used  in  the  numerical  simulations.  Planes  of 
blade-to-blade  computational  grids  such  as  shown  in  Figure  2 
are  stacked  in  the  spanwise  drection  to  create  the  discretized  ge¬ 
ometry.  Computational  grid  lines  are  clustered  adjacent  to  the 
airfoil  surface  and  endwall  to  properly  resolve  the  viscous  flow. 

Figure  3  illustrates  the  predicted  fine  grid  (Case  6)  and  exper¬ 
imentally  measured  pressure  coefficient  distributions  along  the 
surface  of  the  airfoil  at  the  2.9.  12.5.  25.0.  and  50.0%  span  loca¬ 
tions,  where  the  pressure  coefficient  is  defined  ws: 


and  U i  is  the  total  inlet  velocity.  Experimental  data  was  taken 
on  two  separate  blades  in  the  cascade  and  both  sets  of  data  are 
presented  for  each  spanwise  location.  There  is  close  agreement 
between  the  predicted  results  and  the  experimental  data  for  the 
entire  span  of  the  cascade.  In  both  the  predicted  results  (Cases  4- 
6  in  Table  1)  and  the  experimental  data  the  effects  of  the  endwall 
secondary  flow  on  the  pressure  distribution  are  visible  on  the 
suction  surface  of  the  airfoil  near  the  trailing  edge  at  12.5%  span. 
This  phenomenon  was  not  present  in  the  coarser  grid  calculations 
(Cases  1-3). 

The  predicted  and  experimental  [2]  endwall  static  pressure 
coefficient  contours  are  illustrated  in  Figures  4  and  5.  respec¬ 
tively  Comparison  of  the  two  figures  show  that  the  predicted 
results  and  experimental  data  agree  very  well.  Both  figures  show 
a  pressure  ■‘island",  corresponding  the  minimum  pressure  point, 
located  on  the  suction  surface  at  approximately  60%  axial  chord. 
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Figure  3:  Airfoil  surface  static  pressure  coefficient  distribu¬ 
tions  22 

The  predicted  static  pressure  coefficient  contours  for  the  suc¬ 
tion  and  pressure  surfaces  of  the  airfoil  are  shown  in  Figures  6  and 
7,  while  the  experimental  contours  for  both  surfaces  are  shown 
in  Figures  8  and  9,  respectively.  It  is  important  to  note  that  Fig¬ 
ures  6  and  7  are  three-dimensional  renderings  of  the  airfoil  sur¬ 
face,  while  Figures  8  and  9  are  two-dimensional  projections  of  the 
airfoil  surface.  There  is  good  agreement  between  the  predicted 
results  and  experimental  data  on  both  surfaces  of  the  airfoil. 
The  static  pressure  island  associated  with  the  minimum  pres¬ 
sure  point,  located  at  38%  axial  chord  on  the  suction  surface, 
is  shown  in  both  the  predicted  and  experimental  results.  The 
effects  of  the  endwall  boundary  layer  are  visible  in  Figures  6  and 
8  just  downstream  of  the  mid-axial  chord  location.  In  both  the 
predicted  results  and  the  experimental  data,  the  static  pressure 
changes  very  little  for  the  first  50%  axial  chord  of  the  pressure 
surface  (see  Figures  7  and  9),  then  decreases  very  rapidly  as  the 
flow  accelerates  sear  the  trailing  edge. 

The  predicted  limiting  streamlines  for  the  suction  surface  of 
the  airfoil  are  shown  in  Figure  10.  Limiting  streamlines  obtained 
experimentally  using  surface  Sow  visualization  techniques  are 
included  for  companion  in  Figure  11.  The  predicted  limiting 
streamlines  (Figure  10)  show  close  agreement  with  the  experi¬ 
mental  data.  The  flow  in  the  mid-span  region  of  the  suction 
surface  is  two-dimensional,  while  the  flow  in  the  endwall  region 
is  highly  three-dimensional.  A  separation  streamline  forms  at  the 
airfoil  leading  edge  in  the  endwall  region  and  extends  span  wise 
to  approximately  the  25%  span  position  at  the  trailing  edge.  As 
the  flow  moves  towards  the  exit  of  the  blade  passage,  it  is  forced 
towards  mid-span  by  the  endwall  secondary  flow.  The  predicted 
pressure  surface  limiting  streamlines  are  shown  in  Figure  12.  The 
corresponding  experimental  pressure  surface  limiting  streamlines 
are  presented  in  Figure  13.  The  leading  edge  stagnation  line  is 
well  illustrated  in  Figure  12  and  it  is  evident  from  both  the  pre¬ 
dicted  and  experimental  hmitiug  streamlines  that  the  flow  on  the 
pressure  surface  is  predominantly  two-dimensional  in  nature. 


Figure  4:  Predicted  3-D  endwall  static  pressure  coefficient  con¬ 
tours 


Figure  5:  Experimental  3-D  endwall  static  pressure  coefficient 
contours 


The  predicted  endwall  limiting  streamlines  from  the  fine  grid 
simulation  are  shown  in  Figure  14,  while  the  experimental  end- 
wall  limiting  streamlines  are  shown  in  Figure  15.  In  both  Fig¬ 
ures  14  and  15,  a  saddle  point  is  evident  near  the  pressure  surface 
leading  edge  of  the  airfoil.  In  the  numerical  simulation,  the  sad¬ 
dle  point  location  is  predicted  somewhat  further  upstream  than 
that  observed  in  the  experiment.  This  may  indicate  that  the 
predicted  endwall  boundary  layer  is  slightly  too  thick.  In  Ref¬ 
erence  [2]  it  was  demonstrated  that  the  location  of  the  saddle 
point  moves  upstream  as  the  endwall  boundary  layer  is  thick¬ 
ened.  In  both  Figures  14  and  15,  the  flow  in  the  leading  edge 
region  moves  onto  the  suction  side  of  the  passage  and  intersects 
the  suction  surface  near  the  minimum  pressure  point.  The  flow 
on  the  pressure  side  of  the  passage  migrates  towards  the  suction 
surface,  intersecting  it  in  the  region  downstream  of  the  minimum 
pressure  point. 
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Figure  6:  Predicted  suction  surface  static  pressure  coefficient  con¬ 
tours 
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Figure  7:  Predicted  pressure  surface  static  pressure  coefficient 
contours 


Figure  Experimental  suction  surface  static  pressure  coefficient 
contours 


Figure  9:  Experimental  pressure  surface  static  pressure  coeffi¬ 
cient  contours 


HUB  EUDIALL 


Figure  10:  Predicted  suction  surface  limiting  streamlines 


Figure  11;  Experimental  suction  surface  limiting  streamlines 
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Figure  12:  Predicted  pressure  surface  limiting  streamlines  Fi«ure  13:  Experimental  pressure  surface  experimental  limiting 

streamlines 


Figure  14:  Predicted  endwall  limiting  streamlines 


Figure  15:  Experimental  endwall  limiting  streamlines 


To  assess  the  ability  of  the  present  Navier-Stokes  analysis  to 
predict  heat  transfer,  the  Stanton  number  was  calculated  along 
the  surface  of  the  airfoil  and  compared  with  the  experimental 
data  [2],  The  Staatoa  number  is  defined  as: 

St  -  (£)/(lUPr{h„  -  M)  <«) 

and  n  is  the  normal  distance  from  the  wall,  h**  is  the  astiahatk 
enthalpy  at  the  wall,  and  is  the  enthalpy  at  the  wall.  Figure  It 
compares  the  predicted  suction  surface  mid-span  Stanton  number 
distributions  obtained  with  and  without  transition,  the  predicted 
results  from  the  two-dimensional  simulation,  and  the  experimen¬ 
tal  data.  The  fully  turbulent  calculation  overpredicta  the  Stanton 
number  in  the  first  2 &%  of  the  axial  chord,  demonstrating  the  im¬ 
portance  of  modelling  transition.  As  discussed  earlier,  mid-span 


Figure  16:  Mid-span  suction  surface  Stanton  number  distribu¬ 
tions 
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Figure  IT:  Mid-span  pressure  surface  Stanton  number  distribu¬ 
tions 


transition  was  prescribed,  according  to  the  experimental  data,  to 
be  instantaneous  at  the  25%  axial  chord  position.  The  effects 
of  transition  are  evident  from  the  rapid  increase  in  the  Stanton 
number  at  the  25%  axial  chord  location.  The  Stanton  number  is 
predicted  very  well  ahead  of  transition  but  is  slightly  underpre¬ 
dicted  aft  of  the  transition  point.  It  has  been  found,  based  upon 
a  number  of  two-dimensional  simulations,  that  this  small  dis¬ 
crepancy  is  due  to  the  instantaneous  transition  model  currently 
used  in  the  three-dimensional  procedure.  Very  good  agreement 
was  found  to  exist  between  the  predicted  Stanton  number  dis¬ 
tribution  and  the  experimental  data  in  the  two-dimensional  sim¬ 
ulations  when  transition  was  modelled  over  a  region  extending 
between  25  -  60%  axial  chord  (see  Figure  16). 


A  comparison  of  the  predicted  three-dimensional  mid-span 
pressure  surface  Stanton  number  distribution  with  the  experi¬ 
mental  data  is  presented  in  Figure  17.  The  three-dimensional 
predicted  results  demonstrate  only  fair  agreement  with  the  ex¬ 
perimental  data.  A  series  of  two-dimensional  simulations  were 
performed  to  determine  if  the  differences  between  the  predicted 
Stanton  number  and  the  experimental  data  near  the  25%  axial 
chord  location  was  due  to  transition,  similar  to  the  suction  sur¬ 
face.  The  results  from  these  two-dimensional  simulations  showed 
that  transition  had  little  effect  on  the  pressure  side  Stanton  num¬ 
ber  distribution  in  this  region  of  the  Sow.  Further  invest:*alion 
is  required  to  determine  the  contributing  factor  of  this  difference. 

The  predicted  three-dimensional  fine  grid  suction  surface  Stan¬ 
ton  number  contours  are  shown  in  Figure  18.  The  experimen¬ 
tal  suction  surface  Stanton  number  contours  are  illustrated  in 
Figure  19.  The  predicted  results  demonstrate  good  quantita¬ 
tive  agreement  with  the  experimental  data,  except  along  the  hub 
endwall  downstream  of  the  separation  streamline,  where  the  ex¬ 
perimental  data  indicates  increased  amounts  of  heat  transfer.  In 
Figure  18,  the  variation  of  the  transition  location  with  span  is 
apparent  in  the  leading  edge  endwall  regions.  Both  the  predicted 
results  and  the  experimental  data  indicate  enhanced  heat  transfer 
at  the  airfoil  leading  edge  in  the  endwall  regions  and  downstream 
of  the  endwall  separation  line. 

Figures  20  and  21  illustrate  the  predicted  fine  grid  and  ex¬ 
perimental  Stanton  number  contours  for  the  pressure  surface. 
The  predicted  Stanton  number  distribution  shows  close  agree¬ 
ment  with  the  experimental  data.  Both  the  predicted  results 
and  the  experimental  data  reveal  much  less  heat  transfer  taking 
place  on  the  pressure  surface,  as  indicated  by  smaller  values  of 
the  Stanton  number. 


TRAILING  EDGE  TRAILING  EDGE 


MIDSPAN 


Figure  18:  Predicted  suction  surface  Stanton  number  contours 


Figure  19:  Experimental  suction  surface  Stanton  number  con¬ 
tours 
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Figure  20:  Predicted  pressure  surface  Stanton  number  contours  Figure  21:  Experimental  pressure  surface  Stanton  number  con¬ 
tours 


Figures  22  and  23  present  the  predicted  fine  grid  and  exper¬ 
imental  Stanton  number  contours  for  the  endwall-  In  both  the 
predicted  results  and  the  experimental  data,  the  lowest  values  of 
the  endwall  Stanton  number  are  located  on  the  pressure  surface, 
beginning  just  downstream  of  the  saddle  point  and  extending 
to  the  trailing  edge.  In  both  figures,  a  region  of  increased  heat 
transfer  is  located  near  the  suction  surface  leading  edge  and  in 
the  vicinity  of  the  minimum  pressure  point  on  the  suction  surface. 


In  the  wake  region,  discrepancies  exist  between  the  predicted  re¬ 
sults  and  the  experimental  data.  The  predicted  results  indicate 
a  small  amount  of  heat  transfer  in  the  wake  region,  while  the  ex¬ 
perimental  data  indicates  larger  amounts  of  beat  transfer.  These 
discrepa  ;es  are  probably  due  to  the  use  of  an  O-grid  around  the 
airfoil,  which  precludes  the  use  of  wake  modelling  in  the  algebraic 
turbulence  model.  Other  than  in  the  wake  region,  the  predicted 
results  show  good  agreement  with  the  experimental  data 
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Figure  24:  Gap-averaged  total  pressure  coefficient  distributions 


experimental  data.  The  profile  loss  at  mid-span  is  predicted  very 
well  by  the  current  procedure,  but  the  peak  loss  region  near  the 
25%  span  location  is  not  being  accurately  predicted.  The  pre¬ 
dicted  loss  in  this  region  associated  with  the  secondary  flow  is 
predicted  to  be  lower  and  more  towards  mid-span.  The  small 
rise  in  the  predicted  loss  distribution  in  the  endwall  region  ( Cases 
5-G)  agrees  well  with  the  experimental  data. 

Table  1  provides  the  area-averaged  total  pressure  loss  for  each 
of  the  curves  shown  in  Figure  24.  These  values  of  area-averaged 
loss  were  computed  using  Cpj  values  at  the  experimental  data 
locations  in  order  to  be  consistent.  Table  1  shows  that  for  the 
finest  grid  solution,  a  2%  error  exists  between  the  predicted  and 
experimental  area-averaged  losses.  The  primary  source  of  this 
error  is  the  underprediction  of  the  secondary  flow  loss  near  the 
25%  span  location  as  shown  in  Figure  24.  Since  an  essentially 
grid-independent  solution  was  achieved  using  the  finest  grid  (see 
Figure  1),  this  error  can  be  mainly  attributed  to  the  algebraic 
turbulence  model  and  the  transition  model. 


The  predicted  spanwise  distributions  of  the  total  pressure  loss 
coefficient  downstream  of  the  trailing  edge  for  all  of  the  cases 
given  in  Table  1  are  compared  with  the  experimental  data  in 
Figure  24.  The  total  pressure  loss  coefficient  it  defined  as: 


2(iV-fn) 

pUS 


(13) 


The  mid-span  profile  loss  is  accurately  predicted  using  the  fine 
grid,  but  overpredicted  using  the  coarser  grids.  The  predicted 
fine  grid  spanwise  loss  distribution  shows  similar  features  as  the 
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Figure  25:  Predicted  total  pressure  loss  coefficient  contours  at 
exit  of  turbine  passage 


Figures  25  and  26  illustrate  the  predicted  and  experimental 
total  pressure  loss  contours  at  the  exit  plane  of  the  turbine  pas¬ 
sage.  The  effects  of  the  passage  vortex  and  endwall  secondary 
flow  on  the  total  pressure  loss  distribution  are  accurately  pre¬ 
dicted  by  the  numerical  analysis.  Both  the  predicted  and  the  ex¬ 
perimental  total  pressure  loss  contours  indicate  the  existence  of 
three  peak  loss  regions.  The  these  peak  loss  regions  may  be  asso¬ 
ciated  with  three  separate  vortices  which  form  a  triangle  bounded 
by  the  airfoil  suction  surface  and  the  hub  endwall.  The  largest 
vortex,  associated  with  the  secondary  flow,  is  located  at  the  apex 
of  the  triangle  [2].  By  examining  the  vorticity  field,  it  was  de¬ 
termined  that  the  secondary  flow  vortex  rotates  in  the  clockwise 
direction  (see  Figure  25).  A  slightly  smaller  vortex,  which  may 
possibly  be  associated  with  the  suction  side  leg  of  the  horseshoe 
vortex,  rotates  in  the  counter-clockwise  direction  and  is  located 
more  towards  the  mid-span  region  than  the  secondary  flow  vor¬ 
tex  [2].  The  smallest  vortex,  which  may  be  associated  with  the 
suction  side  corner  vortex,  also  rotates  in  the  counter-clockwise 
direction  and  is  located  closer  to  the  hub  endwall  than  the  other 
two  vortices. 

The  predicted  spanwise  distribution  of  the  gap-averaged  flow 
angle  (in  the  pitch-wise  direction)  downstream  of  the  trailing  edge 
for  the  cases  given  in  Table  1  are  compared  with  the  experimental 
data  in  Figure  27.  The  shape  of  the  predicted  three-dimensional 
flow  angle  distribution  demonstrates  good  qualitative  agreement 
with  the  experimental  data,  but  the  numerical  values  are  uni¬ 
formly  larger  than  the  experimentally  measured  flow  angles  by 
approximately  1  degree.  As  shown  in  Table  1.  the  exit  flow  an¬ 
gle  changed  very  little  with  increasing  grid  density,  indicating 
that  grid  dependence  is  not  the  cause  of  the  discrepancy  between 
the  predicted  and  experimental  flow  angles.  The  3.5%  difference 
between  the  predicted  area-averaged  exit  flow  angle  and  the  ex 
penmemal  data  could  possibly  to  be  due  to  inadequacies  in  tbe 
modelling  of  transition  and  turbulence. 
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Figure  26:  Experimental  total  pressure  loss  coefficient  contours 

at  exit  of  turbine  passage  Figure  27:  Gap-averaged  flow  angle  distributions 
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Conclusions 

A  three- dimensional,  implicit,  upwind  Xavier- Stokes  procedure 
has  been  used  to  predict  heat  transfer  and  aerodynamic  perfor¬ 
mance  quantities  in  the  Langston  cascade  for  a  thin  inlet  bound¬ 
ary  layer.  A  gnd  resolution  study  has  been  performed  which 
indicates  that  approximately  900.000  (full  spam  computational 
grid  points  are  required  for  a  single  turbomachinery  blade  row  to 
produce  grid-independent  solutions  for  flows  at  Reynolds  num¬ 
bers  typical  of  those  in  gas  turbine  engines.  In  addition,  this 
study  has  shown  that  overall  aerodynamic  turbomachinery  blade 
row  performance  can  be  predicted  within  29?  using  these  fine 
computational  grids  with  the  Baldwin- Loma.\  algebraic  turbu¬ 
lence  model.  Also,  the  importance  of  modelling  transition  from 
laminar  to  turbulent  flow  within  the  blade  passage  for  the  ac¬ 
curate  prediction  of  heat  transfer  is  demonstrated.  Clearly,  fur¬ 
ther  progress  in  reducing  the  computational  time  required  for 
three-dimensional  Navier-Stokes  procedures  is  required  to  enable 
reliable  and  affordable  grid-independent  solutions  for  turboma¬ 
chinery  blade  row  flows.  This  capability  is  required  before  major 
improvements  in  turbulence  modelling  can  be  made  aDd  three- 
dimensional  heat  transfer  and  aerodynamic  performance  quanti¬ 
ties  can  be  routinely  and  accurately  predicted. 

Acknowledgements 

This  work  was  supported  by  the  Naval  Air  Systems  Command 
under  NAVAIR  contract  N001 10-88-0677  from  the  office  of  George 
Derderian  with  Raymond  Shreeve  of  the  Naval  Post  Graduate 
School  acting  as  technical  monitor  and  the  United  Technologies 
Research  Center  under  the  Corporate  Research  Program.  The 
three-dimensional  simulations  performed  during  this  investiga¬ 
tion  were  performed  on  the  NAS  Numerical  Aerodynamic  Simu¬ 
lation  CRAY-2  supercomputer.  The  authors  would  like  to  thank 
Man  Mohan  Rai  and  Linda  Haines  of  the  NASA  Ames  Research 
Center  for  assistance  in  computational  aspects  of  this  investiga¬ 
tion.  The  authors  appreciate  the  helpful  discussions  with  Om 
Sharma  of  Pratt  A  Whitney  and  Mike  Blair  and  Bob  Dring  of 
the  fluted  Technologies  Research  Center  concerning  interpretion 
of  the  experimental  results. 


References 

[1]  Langston.  L.  S..  Nice.  M.  L..  and  Hooper.  R.  M.,  "Three- 
Dimensional  Flow  Within  a  Turbine  Cascade  Passage," 
ASME  Paper  No.  76-GT-50,  March.  1976. 

[2]  Graztani.  R.  A..  Blair.  M.  F.,  Taylor,  J.  R.,  and  Mayle. 
R.  E  .'  An  Experimental  Study  of  Endwall  and  .Airfoil  Sur¬ 
face  Heat  Transfer  in  a  Large  Scale  Turbine  Blade  Cascade." 
Journal  of  Engineering  for  Power,  Voi.  102,  No.  2.  April. 
1980.  pp.257-267. 

[3]  Chan.  D.  C.  and  Sheedy,  K.  P.,  “Turbulent  Flow  Modeling  of 
a  Three-Dimensional  Turbine,"  AIAA  Paper  90-2124.  July. 
1990. 

[4]  Hah.  C.,  “Numerical  Study  of  Three-Dimensional  Flow  and 
Heat  Transfer  Near  the  Endwall  of  a  Turbine  Blade  Row." 
AIAA  Paper  89-1689.  June.  1989. 

[5]  Hah.  C..  "  A  Navier-Stokes  Analysis  of  Three-Dimensional 
Turbulent  Flows  Inside  Turbine  Blade  Rows  at  Design  and 
Off- Design  Conditions."  Journal  of  Engineering  for  Got 
Turbines  and  Power,  Vol.  106.  April,  1984,  pp.  421-429. 

[6|  Moore.  J  and  Moore.  J.  G„  “Performance  Evaluation  of 
Linear  Turbine  Cascades  Using  Three- Dimensional  Viscous 
Flow  Calculations."  ASME  Journal  of  Engineering  for  Gat 
Turbines  and  Power.  October.  1986. 


[7]  Moore.  J.  and  Ransmayr,  A..  "  Flow  in  a  Turbine  Cascade: 
Part  I  -  Losses  and  Leading  Edge  Effects."  ASME  Journal 
of  Engineering  for  Gas  Turbines  and  Power.  Vol  106.  April. 
1984.  pp.  400-408. 

S'  Lee.  D.  and  Knight.  C.  J  .  “  Evaluation  of  an  O-H  Grid 
Formulation  for  Viscous  Cascade  Flows."  AIAA  Paper  "9- 
0207.  January. 

[9j  Baldwin.  B.  S.  and  Lomax.  H..  "Thin-Layer  Approximation 
and  Algebraic  Model  for  Separated  Turbulent  Flows."  AIAA 
Paper  78-257.  January.  1974. 

[lOj  Roe.  P.  L..  “Approximate  Rtemann  Soivers.  Parameter  Vec¬ 
tors.  and  Difference  Schemes."  Journal  of  Computational 
Physics.  Vol.  43.  pp. 357-372.  19*1. 

[11]  Beam.  R.  M  and  Warming.  R.  F..  "An  Implicit  Factored 
Scheme  for  the  Compressible  Navier-Stokes  Equations." 
AIAA  Paper  77-643.  1977 

[12]  Rai.  M.  M..  “Three-Dimensional  Navier-Stokes  Simulations 
of  Turbine  Rotor-Stator  Interaction."  Journal  of  Propulsion 
and  Power.  Vol.  5.  May-June.  1989. 

[13]  Rai.  M.  M.  and  Madavan.  N.  K..  “Multi- Airfoil  Navier- 
Stokes  Simulations  of  Turbine  Rotor-Stator  interaction." 
AIAA  Paper  88-0361.  1988. 

[14]  Rai.  M.  M..  "An  Implicit.  Conservative.  Zonal- Boundary 
Scheme  for  Euler  Equation  Calculations."  AIAA  Paper  65- 
0488.  1985. 

[15]  Hung,  C.  M.  and  Buning,  P.  G..  “Simulation  of  Blunt-Fin 
Induced  Shock  Wave  and  Turbulent  Boundary  Layer  Inter¬ 
action."  AIAA  Paper  84-0457.  January.  1984. 

[16]  Hung.  C.  M.  and  MacCormack.  R.  W..  “Numerical  Solution 
of  Three-Dimensional  Wave  and  Turbulent  Boundary  Layer 
Interaction."  AIAA  Paper  78-161.  January.  1978. 

[17]  Vatsa.  V.  N.  and  Wedan.  B  W..  "Navier-Stokes  Solutions  for 
Transonic  Flow  Over  a  Wing  Mounted  in  a  Wind  Tunnel." 
AIAA  Paper  88-0102.  January.  1988. 

[18]  Dhawan.  S.  and  Narasimha.  R..  "Some  Properties  of  Bound 
ary  Flow  During  Transition  from  Laminar  to  Turbulent  Mo¬ 
tion."  Journal  of  Fluid  Mechanics.  Vol.  3.  1958.  pp  418-436. 

[19]  Johnson.  D.  A.  and  King.  L.  S..  “A  Mathematically  Simple 
Turbulence  Closure  Model  for  Attached  and  Separated  Tur¬ 
bulent  Boundary  Layers."  AIAA  Journal.  Vol.  23.  No.  11. 
November  1985.  p.  1684. 


